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I. INTRODUCTION 


Motivation 

Owing to their favorable performance characteristics, composite 
materials have been gaining wide use in commercial, military and 
spacecraft applications. Typical commercial applications Include 
building materials, sports equipment, and automobile parts such as 
drive shafts, engines, leaf springs and wheels. The main applications 
of composite materials, however, have been in the aerospace and 
defense related projects. Aircraft and spacecraft structures are 
typically weight-sensitive and advanced composites with high strength- 
to-weight and sti f f ness-to-wei ght ratios have become almost 
indispensable in such structures. 

Spacecraft structures and equipment mounted in these structures 
are required to operate in hostile environments and to be subjected to 
a wide range of dynamic loads. When structural resonances are 
excited, the dynamic loads can produce excessive vibration in the 
structures and equipment. There is a need for the vibration of the 
structures to decay as quickly as possible so that it will not affect 
the normal working and utility of the structures and missions. 
Vibration can be significantly reduced by increasing the damping in 
the dominant modes through the application of active and passive 
damping control devices. 



Active control can be used to suppress the vibration by the use 
of certain active elements such as actuators which alter the dynamic 
response of the structure. It requires the use of special purpose 
hardware, and real-time control -desi gn algorithms applicable to 
individual structural components. Active vibration control technology 
is still in its primitive stages of development and unfortunately, 
feasible solutions are still not available to alleviate the vibration 
problems of real-life complex space structures. 

In passive control, energy dissipation can be achieved by 
external add-on damping devices such as isolators or constrained 
viscoelastic layers. Yet another source of passive damping in most 
structures is provided by their joints and supports. The use of 
passive damping treatments provides a simple and reliable solution for 
vibration suppression, in limited frequency ranges. A balanced 
combination of both active and passive damping Is, however,' likely to 
be a more realistic and practical solution in future aerospace 
structures than the use of active or passive control alone. 

In any case, a detailed understanding of the Inherent damping in 
the material and the structural system is important. Little 
information is available on material damping, especially for advanced 
composite materials. It is difficult to obtain a good understanding 
of the damping phenomena in a material or structure as there are many 
mechanisms involved. All these mechanisms eventually involve the 
conversion of mechanical energy into thermal energy. For composite 
materials, the damping mechanisms, of course, are more complex because 
several factors will influence the damping capacity of the material. 
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v These Include, but are not limited to: 1) Internal friction or 

hysteresis within each of the constituent materials, 2) Interfacial 
slip at the fiber-matrix interfaces, and 3) fiber orientation and 
s/ length. Hence damping cannot be estimated easily, but instead must be 

measured. 

The influences of common environmental factors such as moisture 
v and temperature have been found significantly to degrade the matrix 

dominated mechanical properties of many polymer composites. When 
moisture is absorbed, the composite gets heavier and the resin matrix 
w expands. Matrix swelling and rapid heating may eventually lead to 

surface cracking. The amount and type of degradation that this can 
cause in the mechanical properties of the composite depend on matrix 
'v material, fiber content and orientation, cure quality, service 

temperature, humidity, duration of exposure and geometry of the 
component. Hence, it Is extremely difficult to make theoretical 
predictions regarding the effects of the above environmental factors 
on the mechanical properties, particularly the damping capacity of 
composite materials. 

v Joints are used in deployable structures as well as in structures 

which must be assembled in space. Unlike metals, fiber composites 
cannot be welded together. The alternatives Include the use of 
mechanical fasteners, such as bolts and rivets, or adhesive bonding. 

< These techniques produce joints with vastly different properties. 

Adhesive bonding is preferred over conventional fastening techniques 
w due to a greater uniformity in load distribution, elimination of drill 

holes, reduced weight and increased processing ease. 


It is believed that bonded (and bolted) joints act to enhance the 
damping capacity of the structure. The energy dissipation In a joint, 
however, is a complicated process Involving several mechanisms. The 
characterization of a joint presents many experimental and analytical 
challenges. Because of the complexity of the problem it appears that 
very little work is being done in evaluating the damping in joints and 
built-up structures. Hence there is a considerable need for the 
contribution of knowledge In the understanding of the damping capacity 
of joints. 

Scope of the Study 

The Hubble Space Telescope, which will be placed in orbit by NASA 
in 1989, is a new telescope which will be able to peer far out in 
space and back into time, producing Images of unprecedented clarity of 
galaxies, planets and stars billions of light years away from the 
earth. The Space Telescope requires precise pointing accuracy and its 
optical system is held together by a truss made of graphite epoxy 
composite material. Graphite epoxy is a strong, lightweight composite 
material, developed by NASA in the early 1970s, that expands and 
contracts very little In response to temperature changes. The 
metering truss is 5.3 m long and 2.9 m in diameter and has a mass of 
only 114 kg [1] . 

The pointing control system of the telescope consists mainly of 
optical sensors, a digital computer, gyros and four massive reaction 
wheels. Pointing torques are applied with the reaction wheel 
assemblies. The reaction wheels produce vibration disturbances when 



they are rotated due to Imbalance of the rotor and Imperfections In 
the spin bearings. The vibration disturbances could excite the system 
resonances which would lead to significant truss system displacement 
amplitudes. These would result In the blurring of the optical system 
image of the space telescope. Vibration suppression of the truss 
system is provided by passive isolation systems in the form of viscous 
dashpots [2] . A specific knowledge of the material and system damping 
of the truss system is mandatory in the design of the passive 
isolation systems. Hence one of the main objectives of the present 
research was to evaluate systematically, the damping capacity of the 
graphite epoxy material of the telescope truss system. Furthermore, 
since the telescope will operate In space, it is necessary to 
determine the damping capacity of these materials In a simulated space 
environment. 

When graphite epoxy materials are used In a vacuum environment, 
they have a tendency to outgas. During outgassing they release all 
the gases and moisture which were trapped during manufacture. There 
is some concern, however, that mechanical properties Including its 
material damping may change after the material has outgassed for long 
periods of time. For graphite epoxy materials, this outgassing can 
take a period of about six months but this period may be reduced to 
about 3-4 weeks by heating the material to about 100°C in a vacuum 

chamber. 

The telescope is expected to experience variations in temperature 
of as much as 140°C in space[3]. As mentioned previously, the damping 
capacity of graphite epoxy materials is considerably influenced by 


the variation in temperature. Though it is known that temperature 
affects the damping capacity of these materials, not much work has 
been done on the behavior of these materials in cryogenic 
temperatures. Such temperatures may be important during launching and 
deployment in space. 

The truss system of the telescope is a joint-dominated structure. 
As discussed earlier, the damping provided by joints is a potential 
source of passive damping which is useful for the suppression of 
vibration in the truss system. It is not known, however, whether the 
bonded joints in the truss will significantly increase the damping 
over the material damping. Analytical and experimental work on the 
evaluation of system damping of basic structural joints is a good 
starting point for research in this area. The data derived from 
simple structures will provide some necessary information which can be 
used to gain an understanding of the dynamic behavior of complex 
structures such as the truss system of the telescope. 

The other objectives of the present study include, i) a study of 
the effects of outgassing (moisture desorption), and temperature (both 
high and low) on the damping capacity of the graphite epoxy composite 
material of the truss system, and b) an analytical and experimental 
evaluation of the damping capacity of simple adhesively bonded 
structural joints. 

Report Outline 

This report contains the results of the present 
investigation into the above-mentioned objectives. The project work 


was performed In the Sound and Vibration Laboratory of the Mechanical 
Engineering Department at Auburn University. The report is divided 
into seven major chapters. A general overview of damping— its 
characterization, and some mathematical models Is presented in chapter 
2. The various experimental and estimation techniques for the 
evaluation of damping are included in the same chapter. The advantages 
and disadvantages of each method are discussed. An improved forced 
vibration method developed during the course of the project is 
explained. 

Chapter 3 covers the previous pertinent work on the evaluation of 
damping of composite materials. Work on the effects of moisture 
absorption and desorption, and temperature on the dynamic properties 
of composites is also reviewed. Recent developments in the study of 
dynamics of structural joints is also discussed. 

Chapter 4 deals with the measurement of the damping of graphite 
epoxy material. It details the methodology and results of the damping 
measurements made on graphite epoxy tube and beam specimens under 
normal atmospheric conditions and in a vacuum. Experimental findings 
regarding the influence of moisture desorption and temperature (both 
high and low) on the damping capacity of the material are also 
discussed. 

Chapter 5 is devoted to the study of damping of bonded structural 
joints. An analytical model has been developed to predict the natural 
frequencies, loss factors, and mode shapes of a bonded lap joint for 
free flexural vibration. The mathematical details, solution scheme, 
and numerical results are presented. Numerical results are compared 


with the experimental results obtained on graphite epoxy lap-jointed 
specimens. Some experimental results on the damping of 
graphite epoxy double-butt-jointed specimens are described in the same 
chapter. A summary of the project results, conclusions drawn, and some 
recommendations for future work are presented in chapters 6, 7 and 8. 



II. DAMPING OVERVIEW 


The term damping refers to any form of energy dissipation from a 
vibrating system. There are two main kinds of damping: 1) material 
damping and 2) system or structural damping. Material damping is the 
damping Inherent in the material and Is caused by a variety of 
physical mechanisms which depend on the material. System damping is 
due to energy dissipation in the total structure, i.e. in addition to 
the damping present in the material, it includes several non-material 
damping mechanisms such as, energy dissipation effects at joints, 
'fasteners, and Interfaces, acoustic radiation damping, and coulomb or 
dry friction damping, etc. A knowledge of the damping property of a 
structure is of primary Importance In controlling its resonant 
response, and thus in prolonging its service life. Damping also plays 
a crucial role In the proper design of the system for vibratory 
loadings. 

Characterization of Damping 

A number of different units and notations have been proposed to 
express the damping parameter. The variety of combinations of energy 
dissipation mechanisms, the wide range of materials and testing 
techniques, the effects of joints and interfaces, and various 
different motivations for damping studies have led to different 
viewpoints towards damping and different mathematical models for its 


char ac ter iz at ion [4 -10]. Many different disciplines have been 
concerned with the damping research and each discipline, it appears, 
has its own favorite nomenclature. It is therefore essential to 
review at least some of the important damping representations in 
current use and to indicate the relationships between them. 

Absolute Energy Units 

In terms of the absolute energy units, which are applicable to 
both linear and non-linear materials, the damping energy dissipated in 
the specimen in a cycle per unit volume (U dV ) is a measure of the 

material damping. Thus the total damping energy U^, is 

U D - I v U dv dV, (2.1) 

where V is the volume. The usual unit of U Q is joule/cycle and of 

U dV , joule/cm 3 -cycle. 

The absolute energy dissipation is difficult to measure, In 
practice, in most cases. Hence relative energy units are more 
popular. 

Relative Energy Units 

Relative Energy Units, also applicable to both linear and non- 
linear materials are dimensionless ratios of damping energy to strain 
energy. Four different forms of relative units appear in the 
literature: 

1) Th e loss factor r? is defined as the ratio of energy dissipated in 
a cycle to the maximum energy stored during that cycle. For a simple 


harmonic excitation, f(t) = F sin ut, the response can be described as 
x(t) = X sin (wt + $), where <j> is the phase angle between the force 
and the displacement. The energy dissipated D is equal to the work 
done in that cycle: 

D = J 2ff/w f(t) dt (2.2) 

= t X F sin^. 

The strain energy U in the system is nothing but the energy supplied 

to the system, which is given by, 

U = \ J T f(t) x (t) = \ FXcos^ . (2.3) 

1 0 

Therefore, from Equations (2.2) and (2.3), the loss factor is given by 

1 = 2iU = tan ^* (2 * 4) 

2 ) The quality factor Q which is defined as the reciprocal of the 
loss factor t] is also sometimes used as a measure of damping, 

Q = - . (2.5) 

V 

3) A most popular representation of damping Is the Damping Ratio (£■) . 
It is the dimensionless ratio associated with a resonance frequency 
(modal frequency) and is defined by 



where c is the viscous damping coefficient and c c Is the critical 

damping coefficient, for which the system will cease to oscillate when 

t and released. The damping ratio £, often 


displaced from res 


expressed as a percentage gives a quick visualization of the damping 
present in the system. Values of £ <1 suggest that it is an 
underdamped system, £ = 1 represents a critically damped system, 
and £ > 1 represents an over-damped system. Throughout this report, 
the percent critical damping ratio £ is adapted as the measure of 

damping. 

4) Complex Modulus approach - In many polymer and composite 
materials, the Young's modulus(E) can be treated as a complex 
quantity, since there will be a phase difference between the stress 
and strain[10] . In this approach, we define the following relationship 
between stress a and strain e: 

a = Ee + aE || , (2.7) 

where E is the Young's modulus and a is a constant. Assuming a time 
dependance of the form, e lwt , we can write, a = <r 0 e 1wt , and 

e - € q e 1wt . Then Equation (2.7) becomes 

a 0 = E (1 + iwe) e 0 • (2.7a) 

In Equation (2.7a), the term E (1 + i m) can be treated as a complex 
modulus of the material, which is written in a more convenient 

form," 

E*= E' + i E\ (2.7b) 

where E' = E, the real part is called the 'storage modulus', and E" = 
waE, the imaginary part is called the 'loss modulus'. The tangent of 


the phase angle between the stress and strain is defined as loss 
tangent or loss factor. 

tan^ = T} = • (2.8) 

All of the above relative energy unit representations of damping are, 
however, related to one another by: 

7j = tan^ * 2£ = ^ . (2.9) 

It should be noted that Equation (2.9) is valid for the damping 
measured at a resonance frequency and it is assumed that £«1 (lightly 
damped systems) . 


Experimental Methods for ParopInQ Measuren^nt 


There are several methods available for measuring the material 
damping or damping in general. The choice of the technique to be used 
depends on several factors such as whether the material has a small or 
large value of modulus of elasticity, a small or large value of 
damping, the frequency range of interest, temperature, stress state, 
size and shape of the specimen, specimen configuration and excitation, 
available instrumentation, etc. Comprehensive reviews of the basic 
experimental techniques, along with their relative merits and 
demerits, have been published by Plunkett [11], Bert and Clary [12], 
and Chu and Wang [13]. The different methods for the estimation of 
the damping ratio can be broadly classified as: 

1) Time domain methods, 

2) Frequency domain methods, and 

3) Other special methods developed for testing viscoelastic 
materials. 


Time Domain Methods 

As the name implies, the system response data in the time domain 
are utilized in these methods for the estimation of the modal damping 
ratio. There are two such methods: a) Logarithmic decrement method, 
and b) Impulse response method. These are described in the following 


sections. 



Logarithmic Decrement Method 

This is also called the free-vibration decay method. This is one 
of the oldest techniques used for the estimation of the damping ratio. 
In this method, the free vibration decay (displacement amplitude 
history) of a system to a transient input (an impulse) is measured by 
a transducer and recorded by a recorder. From this decay curve, the 
logarithmic decrement 5, which is related to the damping ratio of the 
system is calculated. The theory behind this simple method is briefly 
outlined here. For a single degree of freedom system (Figure 1 (a)) 
the time response of the system due to a transient input can be 

expressed as 

y(t) = A e'^n* sin(w d t + <(>) , (2.10) 

where, A is a constant defined by the initial conditions, £ = damping 
ratio, w n = undamped natural frequency, w d = damped natural frequency, 

l = phase angle, and t = time. Equation (2.10) represents a decaying 
sinusoidal motion (Figure 1 (b)) and provided £ is small, w d » w n and 

if the sine term in Equation (2.10) is unity at some time t 0 , then it 

will again be unity at times, t^ = t Q + 2?rN/w d » t Q + 2 jtN / w n , where N 

is an integer number, N = 1,2,3 Thus, by examining Equation 

(2.10) , we see that, 

N(2lf). (2.n) 

By taking natural logarithms, we obtain the logarithmic decrement 5 


y = disp locemcnt amplitude 


(a) A Single Degree of Freedom Syscem 

(b) Ideal Free Vibration Response when y is Plotted on a 
Linear Scale 

(c) Ideal Free Vibration Response when y is Plotted on a 
Logarithmic Scale 


Figure 1. 


Illustration of the Logarithmic Decrement Method 




17 




V 


w 




/ 


i . 

5 = _ N £n yltJJT 


( 2 . 12 ) 


Use of Equation (2.12) provides a convenient way of measuring the 
damping ratio In this method, the specimen is clamped in a 
cantilever position, then "plucked" and allowed to vibrate freely. 
The decaying motion may be measured by strain gauges, accelerometers 
or optical sensors. Crawley and his co-workers [14,15] have utilized 
strain gages to test their specimen in free-fall in a simulated space 
environment. This was done to avoid the extraneous energy losses due 
to supports and acoustic radiation. 

If the system has many resonances, this approach can still be 
used. One way is to use a band-pass filter to filter out the 
undesired frequency content of the signal prior to displaying the 
signal in the time domain on the recording system. An alternative 
method is to 'tune 1 the desired mode of vibration by attaching several 
shakers to the system and by using force appropriation techniques 
described by Asher [16]. When the system is vibrating in a desired 
mode, the shakers are simultaneously cut-off, and the resulting decay 
for that mode is measured and the damping ratio is calculated by using 
Equation (2.12). This approach, however, requires elaborate 
•instrumentation and hence has been used only for testing large 
aircraft structures. 

In the logarithmic decrement method, it is assumed that the 
damping ratio is independent of the displacement amplitude. A 
modification of the decay method can be used very easily by passing 


the time response through a log amplifier. The decay curve will then 
be a straight line and the slope of this line is a measure of the 
damping [Figure 1(c)]. An alternative approach is to use a least- 
squares curve-fitting procedure for the measured decay curve to obtain 
an accurate estimate of the damping ratio. The theory behind this 
curve-fitting method is given in [17]. A computer simulation was 
carried out in reference [17] in order to compare the decay method 
with the least squares curve-fitting method for the estimation of the 
damping ratio. It was shown that the curve-fitting method is more 
accurate than the decay method for the conditions most frequently 
encountered in practice. The optimum number of data points for which 
the best estimate of the damping ratio is obtained was found to be 
about twenty. 

Some other methods of measurement associated with the decay of 
free vibration of mechanical systems have been noted in the 
literature[4,18] . One such approach Is the use of the time-averaged, 
squared response y 2 , in decibels. If we define the displacement level 
Ly to be: 

Ly = 10 log 10 (y 2 /y 2 ref ). dB (2.13) 

where, y is the displacement amplitude and y^g-f a reference 

displacement amplitude, then it is easy to show, by substituting for y 
from Equation (2.10), that, 


(2.15) 


and the Decay rate A t is given by: 

-dL w 

A t = = - [20 Log 1() (e)] [-£»„] • 8.69^ n . 

Another way of estimating the damping ratio of a system is to 
measure the reverberation time of the vibration decay In the system. 
The reverberation time T R is defined as the time in which the 

displacement level is reduced by 60 dB. It is seen from Equation 

(2.15) that this value is given by: 

T r = 6.90/^ n = 1.10/£f n . (2.16) 

(This is easily shown by letting dt = T R and dL y = -60 dB, and w n = 

2jrf n in Equation (2.15).) 

Impulse Response Method 

This method is essentially a forced vibration method. With the 
advent of modern digital FFT analyzers, it is now possible' to obtain 
the frequency response of a system to any physically realizable 
excitation, and the data can be processed either in the time domain or 
in the frequency domain within a fraction of a second. The frequency 
response H(w) of a system is the ratio of system response to an input 
excitation in the frequency domain. The impulse response h(t) is the 
time domain equivalent of the frequency response. It can be obtained 
very easily by an Inverse Fourier transformation of the frequency 
response function (FRF). 

A typical experimental set-up for measuring the impulse response 
function is shown in Figure 2. For a single degree of freedom system, 
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the unit impulse response (l.e. the system's response to a unit force) 
Is identical to that of Equation (2.10) and it can be easily 
recognized that the exponential decay Is also a straight line when the 
- displacement amplitude Is plotted on a logarithmic scale. Hence, the 

techniques described in the previous section can be used to estimate 
the damping ratio. The decay rate method discussed previously is most 
, commonly used in conjunction with the Impulse response method. 

For a multi-degree of freedom system, the impulse response 
function Is more complicated. It can be written In terms of the modal 
parameters as [19,20]: 

h(t) = [ |R.- |e“^i w i t sin(«.t + a.) , (2.17) 

1 = 1 1 

* where, n = number of modes 

R^= residue, 

< 2 ..= angle of residue, 

£,= modal damping ratio, and 

modal frequency. 

„ The subscript i stands for the ith mode of vibration. The residue R^ 

is related to the mode shape of the ith mode. Brown et al . 1]20] have 
described a curve-fitting technique called the complex-exponential 
y algorithm, to extract the modal damping ratios and other modal 

f parameters from the measured impulse response data using Equation 

(2.17). This technique has the disadvantage that it requires very 
elaborate data reduction equipment and generally gives an 





overestimation of the damping values because of the damping Induced by 
the exponential window in the FFT analyzer. This method is, however, 
the only powerful technique available for measuring the damping ratio 
of systems with very high modal density and modal interference. Many 
commercial modal analysis software packages have Incorporated this 
technique as one of the standard curve-fitting routines for extracting 
the modal parameters. 

Frequency Domain Methods 

Several techniques are available for estimating the damping from 
the system response data in the frequency domain. Usually, in all 
these methods, a simple single degree of freedom system model with 
either a viscous or a hysteretic type of damping is assumed. A 
detailed description of viscous and hysteretic damping models can be 
found in [10]. The following is a broad classification of the 
frequency domain techniques for the estimation of the damping ratio: 

1) Resonant forced vibration methods, 

2) Non-resonant forced vibration methods and, 

3) Input-power method. 

Resonant Forced Vibration Methods 

As in the time domain methods, the damping ratio is assumed to be 
independent of the displacement amplitude, but dependent on the 
resonance frequency of vibration. 
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If the damped simple system in Figure 1 (a) is excited by a 

simple harmonic force of amplitude | F | at angular frequency u - 2/rf, 

the frequency response function is given by[18]: 

Y (w) _ 1 _ 1/K (o 

F(wT ' (K-Mw 2 ) + lew (1-r 2 ) + i(2^r)’ 1 * ' 

Here r = u/w^ is the frequency ratio, and 1 = (-1) ^ 2 . The ratio Y/F 

is called the frequency response function (FRF) which in this case is 
the receptance or the displacement admittance function. In practice, 
transient, random, and sinusoidal excitations can all be used to 
obtain the FRF function. When random excitation is used, the FRF is 
estimated by power spectrum averaging. In this case the FRF is 
defined as the ratio of the cross-spectrum between the force and the 
displacement signals to the auto-spectrum of the force signal. The 
random excitation technique is normally used to obtain more accurate 
results than the other excitation methods if noise Is pres_ent in the 
measurement process. A survey of the different excitation techniques, 
their relative merits and demerits for modal analysis of structures 
can be found in references [21,22]. 

Equation (2.18) forms the basis for many of the damping 
estimation methods as follows: 


(a) Half-power Points Method 

The magnitude of the receptance function of Equation (2.18) can 


be expressed as: 



_ (1/K) j / 

" ((1-r 2 ) 2 + (2£r) J ) 7 2 * 


(2.19) 
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Provided the damping Is small, which Is the situation In most 

engineering applications, then, when the forcing frequency u equals 
the undamped natural frequency w n , the displacement admittance is 

nearly maximum. The dynamic magnification factor, DMF - | Y | / | F | 
/ K for this situation when u = w n (resonance) is thus, DMF = 1/2$- = Q 

(where Q is the quality factor). 

Defining the half-power points as the points at which 


| Y | / | F | / K has fallen to V2 of its value at resonance (see 
Figure 3(a), we obtain the two values of angular frequency Wp and u 2 


for this to occur from 


1 = 1 

2$- V2 { ( 1 - r 2 ) 2 + (2 £t) 2 } 1/j 

and assuming £*«1, we obtain from the quadratic solution: 


( 2 . 20 ) 


Thus, the damping 
bandwidth Aw = Wj-Wi 




2 


LJ 


1 2 
, l . 


i * r* 


( 2 . 21 ) 


ratio £• can be obtained from the half power 
as 


Wj -Wj 

£ = ~2uT~ 


( 2 . 22 ) 


Thus, by determining the two frequencies u 1 and w 2 , at which the 
response is 1/V2 of the value at the resonance frequency w n , the 


damping ratio £ can quickly be determined using Equation (2.22). The 
amplitude ratio of 1/V2 corresponds to reduction of amplitude measured 





In decibels of 20 Log 1Q (l/V2) = -3 dB. Thus a measurement associated 

with an amplitude ratio of 1/V2 is frequently referred to a "3 dB 
bandwidth." If the system has many higher resonance frequencies, then 
the same procedure can still be used by measuring the half-power 
bandwidth of each resonance peak in the curve. This approach can be 
used provided that £ is not too large and that the resonance peaks are 
fairly sharp and can be separated in the frequency domain. 


(b) Damping from Real and Imaginary Parts of Receptance (co-quad 
plots) 

The real and imaginary parts of Equation (2.19) shown in Figure 3 

(c) and 3 (d) are 


I > - ■ and 

( F ) = TT^tW’- 

The real part will be zero at a resonance (when r = 1). 
ratio can be estimated from the relation: 

- u\ 

£ ~ 2 (w 2 2 + u [ 2 ) 1 


(2.23) 

(2.24) 

The damping 


(2.25) 


where u l and w 2 are the frequencies on each side of u n where the 


magnitude of the real part reaches a maximum. 

From the imaginary plot, the resonance is located at the point 
where Im(Y/F) has a peak. The half-power points correspond to 
frequencies at which the quadrature response (imaginary part) has half 
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of the maximum amplitudes. Thus Equation (2.22) can once again be used 
for calculating the damping ratio. 


(c) Damping Estimation from the Phase Angle Plot 

Estimation of the damping ratio from the phase angle plot Is well 
described by Pendered and Bishop [23]. The equation for calculating 
the phase angle is, 

tan ?> ■ • ( 2 - 26 ) 


A typical variation of the phase angle with frequency is shown In 
Figure 3 (b) . At a resonance, theoretically, the phase angle between 

force and displacement should be -90°. The half-power points, A and B 

correspond to phase angles of -45° and -135 c respectively, thus, the 
damping ratio can be evaluated using Equation (2.22). In practice, 
since it Is often very difficult to measure the half power points from 
the phase angle curve, an alternate approach has been suggested. This 
Involves a measure of the slope of the tangent to the curve at the 
inflection point, a 2 ^/9w J = 0. The damping ratio can then be 
calculated from the relation: 


1 


1 


V= 7“ T = (2.27) 

^ n (slope at «„) «„( 5 J ) 

u “n 

Although the phase method has an advantage in the sense that, the 
determination of the natural frequency does not depend on the accurate 
location of a 'peak', it has not, however, been very popular, since, 
even now in practice, it is very difficult accurately to measure the 
phase angles. 




(d) Vector Diagram Method (Kennedy-Pancu Method) 

This method was originally developed by Kennedy and Pancu [24]. 
If the complex receptance function of Equation (2.18) is plotted on an 
Argand plane, which is a plot of real part vs. imaginary part on a x-y 
coordinate system, the resulting curve would be a circle as shown In 
Figure 4. In fact, this can be proved very easily by considering 
Equations (2.23) and (2.24). By a simple mathematical manipulation, 
it can be shown that 

(Re(Y/F)} 2 + {Im(Y/F) + l/(4£r)} a = l/(4^r) a , (2.28) 

which represents a circle with center at (0, 1/(2£t)), and having a 
radius of l/(2£r). From the plot, Kennedy and Pancu noted that: 1) at 
resonance, the displacement vector lies on the imaginary axis and 
hence, Is 90 degrees outof phase with the applied force, 2) the 
change. in arc length(s) along the curve per unit change In frequency 
(ds/dw) Is a maximum at resonance, and 3) the diameter of -the circle 
is Inversely proportional to the damping ratio. 

As in the phase angle method, here also two approaches are used 
to estimate the damping ratio. One approach, of course, Is to use the 
half-power point Equation (2.22). The half-power point frequencies 
can be determined readily from the Argand Plot by drawing a line 
parallel to the real axis through the center of the circle; the points 
of intersection of this line with the circle are the half power point 
frequencies. These points also correspond to frequencies w 1 , and at 

±90° angle with respect to the damped natural frequency. The 
resonance frequency is the point on the circle farthest from the real 


axis (where the Imaginary part is zero). An alternate approach to 
determine the resonance frequency Is from arc length measurements; the 
resonance frequency corresponds to a point of maximum rate of change 
of arc length with frequency. A mathematical analysis for this is 
presented in many references [25-28]. Alternatively, the damping 
ratio, can be estimated by knowing the rate of change of phase with 
frequency, evaluated at resonance (see Equation 2.27). 

In applying the Kennedy-Pancu method, the normal procedure Is to 
fit the 'best circle' through the measured data points around a 
resonance. This is called the 'circle fit' procedure. For a multi- 
degree of freedom system, the complex receptance can be viewed as a 
summation of the complex receptances of several single degree of 
freedom systems as in Equation (2.18). A typical Kennedy-Pancu plot 
for a two-degree of freedom system is shown in Figure 4 (b) . It is 
evident that there are two resonance frequencies because of the 
presence of two circular arcs. The normal procedure for such systems 
is to fit a circle to each of the loops separately as a single degree 
of freedom system. Notice that the vector diagram of systems with 
many degrees of freedom are not circles, but curves with many loops, 
usually one for each resonance. Woodcock [29] has extended the 
Kennedy-Pancu method to study systems with many degrees of freedom. 
It is reported in the literature that the circle-fitting procedure 
yields better estimates of modal parameters than the other methods in 
the presence of closely spaced modes. But, the choice of data points 
utilized in the circle fit gives different answers and the best answer 
becomes a matter of judgement [20]. 



The principal assumptions made in all the methods discussed so 
far for the estimation of the damping ratio are: 1) the system is 
linear, 2) the system possesses very light damping, and 3) the modes 
of the system are well separated. Each method has Its own advantages 
and disadvantages. In general, most of the above methods have a major 
drawback in that they just use two or three points on the response 
curve for the estimation of the damping ratio. Hence, the accui acy of 
the estimated damping ratio depends on how well the two or three data 
points are chosen. In practice, it is usually difficult accurately to 
locate the resonance frequency and the half-power points. 
Furthermore, the frequency resolution of the analyzing Instrument has 
a restriction on the value of the damping ratio that can be estimated. 
Hence, recognizing the need for a better method for damping ratio 
estimation, the following improved forced vibration method has been 
developed as a part of this research project. 

Improved Resonant Forced-Vibration Method ? 

The method consists of fitting a ’best curve' for the measured 
receptance data using an iterative least-squares error criterion. The 
damping ratio and the undamped natural frequency are then, computed 
from the coefficients of the rational fraction polynomial determined 
from the curve-fitting technique. Any of the three functions of the 
receptance namely, the magnitude, the real part, or the imaginary part 
can be used in the method. The theory behind this method is as 
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The frequency response function used In the analysis Is the 
receptance function (displacement/force) of a single degree of freedom 
system, obtained from a single point excitation. Either a random 
input excitation (white noise) or a swept sine excitation can be used. 
The receptance function can easily be measured in practice by an 
impedance head and a dual channel signal analyzer. In the following 
analysis, the equations of motion of a multi degree of freedom system 
with proportional damping are used. These are chosen to illustrate the 
fact that the method developed here can also be used for a system 
response containing several distinct modes by fitting each of the 
modes separately as a single degree of freedom system. 

For a general multi-degree of freedom system with proportional 
damping the equation of motion in matrix form is. 

Mx + RX + Kx = P(t) , ( 2 . 29 ) 

where, M ■ mass matrix, nxn, 

R = damping matrix, nxn, 

K = stiffness matrix, nxn, 

x , X, x = acceleration, velocity and displacement vector, 
respectively, nxl, 

P(t) = force vector, nxl, 
n = number of degrees of freedom. 

If the normal mode shape of the system without damping is [^] nxn » ther 

applying the following transforms: 

x = [^] q f * = [^]Qi anc * x = 
and the orthogonality conditions: 
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[^]j T M[^3-j = 0, [^] i T M[^] i - Mg. 

[*]/«*], = 0, [^] i T K[^] i - Kg, 


where M and K r represent the decoupled (diagonal) mass and stiffness 
matrices of the system respectively. Pre-multiplying both sides of 

Equation (2.29) with [^] T and applying the above transforms the 
following equation is obtained: 

[0] T M[^]q + [^] T R[^]<5 + [^] T K[^]q = [£j T P(t) . (2.30) 

Assuming that the damping is viscous, and that the damping matrix Is 
proportional to mass matrix, r f ■2£ , y .w nr M r , a diagonal matrix equation 

from Equation (2.30) is obtained: 


M G q + Mg[ V 2^ nr \]q + K G q = Pg(t) 


(2.31) 


where, Pg(t) = [^] T P(t). If Pg(t)= Pge 1wt . then the steady state 
solution of Equation (2.31) can be obtained as: , 

q = Qe , 


Gr 


where, Q r = M Gr [wj; r + i2£ r w nr w - u c ] 


2t t r « li 2, i • • n • 


Since x = ty]q = q Q e 


iwt 


n {» r } T M {^r) 

r=1 M Gr[“nr + 1 Vnr“ " A 


(2.32) 
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If only one force acts on the j-th point and the response Is measured 
at the same point, then the r-order driving point displacement 
admittance is obtained as follows: 

,2 


M 




Dr 


(2.33) 


M GrKr + 1 2 £ w nr w " w J 
The magnitude, the real part and the imaginary part of the receptance 
for the r-th mode are, respectively, 




tl /M 4 


+ - i>4f 2 + 


Re(M Dr ) - 


tl (f 2 - f 2 )/(2f) 2 


« Gr [f 4 ♦ 2fe 2 - D^f 2 ♦ f 4 ] ' 


I”( m dP * tL ,,4 


-2# 2 f f„ r f /(2ir) 2 


Gr[f’ r + 2(2f 2 - !)f 2 r f 2 + f 4 ]' 


(2.34) 


(2.35) 


(2.36) 


where, f is the excitation frequency, Hz. Equations (2.34) through 
1 nr 

(2.36) can be expressed by rational functions: 

2 ?3 

l*W = [A x + A 2 f 2 + f 4 ] ' 


(2.37) 


Re(M Dr ) = 


Im(M D ) = 


A 4 (A x ) 1/2 -A 4 f 2 
A x + A 2 f 2 + f 4 

A 5 f 

A 1 + A 2 f2 + f4 


(2.38) 


(2.39) 


where, Aj = f 4 r j A 2 = 2 (2£ 2 -l) f 2 r , A 3 = ^ 2 / M Gr (2jr) 4 , ( 2 . 40 ) 
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A 4 = & M Gr (2^) 2 , and A 5 = 2^/ M Gr (2;r) 2 . 

From the coefficients Aj and A 2 , the natural frequency and the damping 

1 

ratio can be obtained: 


f = A 1/4 
T n fl l ' 

(2.41) 

C • { 1 \\ + m 17 ? 

2 2f 2 

(2.42) 


Using the measured data for the magnitude or the real part or the 
Imaginary part of the receptance and the frequencies in the range of 
the bandwidth of the half-power points, these data can be fitted to 
the corresponding equations given above. Using the magnitude with 
Equation (2.37) for example, the curve-fitting procedure is as 
follows: 

Let y = I M Dr (f ) I ' then the est ’ imate 0 f y 1s: 


^ . B 1 . 
y 1 + B 2 f 2 + B 3 $ 

where, Bj = ^ , B 2 = A^ ' anc ^ B 3 = TTj* ^ c l ua ' t ' Ion (2*43) 
written as, B, - y (1 + B 2 f 2 + B 3 f 4 ) = 0 . The residual r is: 


(2.43) 
can be 


r= J_y=B 1 - Y(B 2 f 2 + B 3 f 4 ) - y (2.44) 

Y = r+y. For each experimental point, corresponding to a y^ - f ^ 


pair, the residual is 


r, ■ Bj - Y, CB 2 ff + B 3 fJ ) -y r 


(2.45) 



The least-squares curve-fitting criterion requires that the sum S of 
the squares of residuals for N points be a minimum, that is 


S * [ r? = [ [B,-Yj (Boff + B 3 f7 ) - y *r s min. (2.46) 
1=1 1 i=l 1 


Setting the partial derivatives of Equation (2.46) with respect to the 
coefficient B^ B 2 and B 3 to be zero, the matrix equation is obtained: 



N 

- E V? 

-EY.f? 

ii 


V 




E Vi 

-EY E f* 



B 2 

= 



JVi 4 

-n? f* 

-EY E f* _ 


__ B 3_ 



Equation (2.47) 

can be solved 

for the constants. 

B l' 

B 2 and B 3 


(2.47) 


by an 


Iterative method. For the first Iteration, It Is assumed that 
Y.j=y.j ,then, (L) = y^ + r^(L-l), where L Is the Iteration time. Then 

the calculations are repeated till the desired convergence driterion 
Is obtained. After at most 10 Iterations, the coefficients will 

converge to an accuracy of about 10 -6 . Then the undamped natural 
frequency and the damping ratio can then be obtained from Equations 
(2.41) and (2.42) respectively. A computer program has been written 
to carry out the iteration using all the three functions, namely the 
magnitude, the real part and the imaginary part of the receptance 
data. 

The residual for Equation (2.43) can be defined in an 


alternative form: 
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r 1 = y i (A l + A 2 f2 + f4) - A 3' (2 ' 48) 

then the matrix equation for the least squares curve-fitting will be 
simpler than Equation (2.47) and this procedure does not require any 
iteration for its solution. Hence this is even a simpler method than 
the Iterative least square method suggested above. Similar expressions 
can be developed using the real and imaginary parts of the receptance. 

The major advantage of this method is that, the entire data near 
a resonance are used for the estimation of the damping ratio, as 
opposed to three points in the half-power points method. Thus, a more 
accurate estimate of the damping ratio can be obtained. Furthermore, 
the undamped natural frequency obtained from the least-squares 
computation is used in the calculation of the damping ratio. Also, 
this method is not limited to systems with small damping and the 
frequency resolution of the analyzer does not Influence the accuracy 
of the method to some extent. Since this method uses fewer data near 
a resonance as compared to other complex curve-fitting techniques 
[30,31] a smaller memory is required in the computer and a higher 
data-reduction speed can be achieved. This method has been used in 
the present investigation for the measurement of the damping ratio of 
the graphite epoxy composite material. Details of its application to 
the present project are discussed in the forthcoming chapters. Some 
results showing the comparison between the curve-fitting technique and 
the half-power points method is Included in Appendix A. 
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Non-Resonant Forced Vibration Method s 

This method has been mostly applied In the measurement of the 
complex Young's modulus of polymers and viscoelastic materials. The 
loss factor ( 2 ^-) is calculated from the real and imaginary part of the 
modulus by use of the equation 7 = E"/E', where E' and E" are the real 
and imaginary parts of the complex dynamic modulus respectively. The 
method can be best illustrated by considering the equation of motion 
for a single degree of freedom system using the complex modulus 
approach. For many composite materials, the force and displacement 
cycles will be out-of-phase by <f>, so if F Q and y Q are the force and 

displacement amplitudes respectively, then 

M y + (K 1 + iK") y = F Q e 1 + ^ , (2.49) 


where u is the angular frequency. ■ 

Let y = y 0 e iwt , then Equation (2.49) can be written as: 

-Mw ! + K' = (F 0 /y 0 ) cos<p, and (2.50) 

K" = (F 0 /y 0 ) sin^ . (2.51) 


For a very small sample, operating at low frequencies (stiffness- 
controlled region), the inertia term of Equati on (2 . 50) is usually 
neglected. Further, the stiffness K of the beam specimen is related 
to the modulus E in the form K = EA/t, where, A is the area of cross 
section and t the thickness of the sample. Hence, from Equations 
(2.50) and (2.51), the damping ratio can be evaluated as: 

£ = K"/K' = 1/2 tan <f>. 


(2.52) 


Hence the forced non-resonance method relies on the measurement of the 
phase angle between the force and displacement signals at a particular 
frequency w, thus, It can be used as an Indirect method for obtaining 
the damping as a function of frequency. This method, however, has a 
major drawback in that it can be used only on very small specimens and 
In the low frequency region, since the effect of inertia is neglected 
in the method. This method has been adapted for the measurement of 
the complex modulus of polymers and composites in many commercially 
available Instruments like the DMTA (Dynamic Mechanical Thermal 
Analyzer) [32,33]. 

In a typical DMTA test arrangement, the displacement Is applied 
via a loading stage attached to the table of an electromagnetic 
vibrator. The displacement is measured usually by a non-contacting 
magnetic transducer. The force is measured by a force transducer 
attached behind the specimen. The specimen is mounted horizontally 
within a controlled temperature enclosure in a single cantilever or 
dual -cantilever. A typical sample used in this arrangement is about 1 
mm thick, 6 mm wide, and 28 mm long, and the frequency of excitation 
Is limited to below 90 Hz. 

The advantage of this method lies in its ability to study the 
dependence of the dynamic properties upon frequency, strain amplitude 
and temperature. Materials with a modulus greater than 10 GN/m 2 cannot 
easily be tested by this method. For the measurement of the damping 
the phase angle must be measured accurately. For phase angles less 
than 1 degree it is not possible to measure them very accurately. But 
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larger values of phase angles can be measured accurately. The method 
therefore Is not good for measuring small values of damping, but Is 
Ideally suited for measuring large damping values. 

Yhe hysteresis loop approach is another widely used non-resonance 
forced vibration method. Its theoretical basis is explained in detail 
by Lazan [2]. This method enables the measurement of very low 
material damping particularly in metals. It is limited, however, to 
the very low loading frequencies that are usually achievable in 
material testing machines. 


Input Power Method 

The theory for this method is developed from a consideration of 
the energy balance in the structure. If a- damped mass-spring single 
degree-of-freedom system (Figure 1 (a) ), is excited by a force F at a 
sinusoidal frequency w, then the dissipated power n d .j SS 

proportional to the product of damping force F d and velocity: 


n 


di ss/cycle 


2?r/w 

/ F. dv = / F. v dt = ncuv ( . 


(2.53) 


where v is the velocity amplitude, and c is the viscous damping 
o 

coefficient. 

In steady state conditions, the power supplied II^ n from a shaker 

must be equal to the power dissipated. Thus If can be measured 

and II- = IIj- is assumed, then c can easily be calculated from 
in d i s s 

Equation (2.53). The input power can easily be measured using an 


• v 


w 
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Impedance head. If a simple harmonic force of angular frequency u and 
amplitude F Q is applied, then the time-averaged input power is 

n in " v o F o cos ^ (2-54) 

where p is the phase angle between the force and velocity. Fahy [34] 

has discussed the case of input power measurement for random 

excitation. Nakayama and Crocker [35] have made extensive 

measurements of power supplied to structures and determined the 

structural damping from such measurements. According to Nakayama and 

Crocker, the input power for a random excitation with a time delay 

between the force and acceleration signals is given by: 

N/2 f k + Af/2 

n i n ■ L (C Fa (f)sinA^ k + Qp a ( f k ) c °sA^ k } • *n [ f _ Af/2 ],(2.55) 
K= 1 

where, Af ■ frequency resolution of the analyzer, 

N = number of data points in the time domain (Block size). 

C_ (f. ) and Qr (fj are the real and Imaginary parts of the cross 

spectrum between force and acceleration respectively at frequency f^. 

A^ k = phase lag. between acceleration and force signals. The band- 

limited input power can be obtained from Equation (2.55) by a simple 
averaging procedure. Structures such as panels have many resonance 
frequencies, and the normal procedure, then, is to measure the space- 
averaged mean square velocity over the structure: < vJ > Sjt - The 

dissipative loss factor 77 = 2£ is then obtained as [36]: 

^in 

^ = M<v 2 >2?rf 1 


(2.56) 


where M is the mass of the structure. The input power method has been 
used mainly in conjunction with an approach called the Statistical 
Energy Analysis of Structures (SEA) [36] . Cremer, Heckl and Ungar [37] 
have developed an expression for the loss factor based on the modal 
density (number of modes in a frequency band) for a broad band random 
excitation. Clarkson and his co-workers [38-40] have used these 
equations to obtain the modal densities and frequency-averaged loss 
factors of plates and shells. The frequency-averaged loss factor 
results they obtained compared well with the results from the 
logarithmic decrement method. According to their experience, at least 
five natural frequencies should be present within the analysis band 
and the damping should be small for the results to hold good. 

The input power method is only an approximate method capable of 
yielding frequency-averaged loss factors. 

Other Methods 

Several other methods for measuring the damping ratio have been 
developed specifically for viscoelastic materials. Some of the 
Important techniques are briefly reviewed here. The Oberst bar 
technique - also called the Bruel & Kjaer cantilever beam method, was 
developed by Oberst [41] and is the basis for the ASTM standard method 
for measuring the material damping [42]. In the original Oberst 
method, a cantilever metal beam (steel or aluminum) is coated on one 
side with a layer of viscoelastic material whose damping is to be 
measured. The beam is excited by a sinusoidal signal with different 
frequencies and the frequency response of the system is measured. The 
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loss factor t ) and the Young's modulus of the material are then 
extracted from lengthy equations given In [45]. Due to the high 
thermal coefficient of expansion for most viscoelastic materials, 
considerable thermal bending is induced on the coating side of the 
Oberst bar. In order to cancel out the thermal bending, Nashif [43] 
proposed to coat both sides of the beam. He derived separate 
expressions for calculating the modulus and the loss factor. In this 
method, the loss factor of the metal beam is assumed to be negligible. 
Jones [44], further improved the testing apparatus with an approach 
that allowed the test specimen to be heated to any temperature within 
the limitation of the testing apparatus. 

Cannon, Nashif and Jones [45] have proposed another method for 
measuring the damping ratio of very soft viscoelastic materials. In 
this method, a cylindrical specimen of the material with a variable 
mass on top is driven by a shaker. The loss factor can be calculated 
by measuring the amplification factor, which is the ratio of the 
resonant displacement response of the top mass to the shaker input. 
Although, this method requires a simple test setup, it is not very 
accurate for materials with modulus greater than 7 MN/m 2 . 

Ritcher [46], has proposed a rotatinq-beam deflection method for 
measuring the damping ratio. In this method, a tubular shaft, made of 
the material under investigation is mounted in a cantilever mode and 
rotated about its axis in a horizontal plane. The horizontal and 
vertical displacements of the free end of the shaft are measured 
optically, and their ratio is related to the loss factor. The 


rotatl ng-beam deflection method, however, has not been popular since 
it requires very-high-precision optical instruments and is limited to 
very low frequencies (0.001 to 20 Hz). 

Damping has been estimated In the special case of a base-excited 
system by measuring the ratio of tip to base displacement amplitudes 
of a cantilever-beam specimen vibrating in its fundamental mode [47- 
50], Gibson and Yau [51], have derived a more general expression for 
calculating the damping in terms of the displacement amplitude at an 
arbitrary point along a beam vibrating In any mode. In a previous 
paper, Gibson and Plunkett [52], have measured the damping ratio of a 
double-cantilever beam specimen from the energy balance method. Under 
steady state conditions, the material damping is defined in terms of 
the ratio of input energy to the strain energy stored in the systems. 
In their experiments the specimen was driven in Its first or second 
mode, and the resonance frequency, input acceleration and bending 
strain (using a strain gage) were measured. The input energy and the 
strain energy were calculated from the above-measured quantities and 
the specimen geometry. The accuracy of the damping ratio obtained by 
this method depends on the accuracy with which the bending strain is 
measured. 

A method utilizing the gravitational acceleration called the 
gravi ty m ethod has been proposed by Sekiguchi and Asami [53] for 
measuring large damping. The authors have used this method for 
measuring the viscous damping coefficient of an oil damper. Several 
types of wave propagation methods have been reported in recent years 



In connection with the measurement of dynamic modulus of long and 
slender composite specimens or structural elements [12,54]. The 
dynamic properties of the test specimen are measured by observing the 
changes that occur in certain characteristics of a stress wave during 
its propagation through that medium. Both pulse and continuous waves 
can be utilized. The main advantage of this method Is its ability to 
generate damping data that are independent of the damping of the 
support system. The drawback, however. Is that tests are confined to 
long and slender rods subjected to high loading rates. 

Summary 

Among the various methods for measuring the damping ratio of 
materials and structures described In this chapter, the logarithmic 
decrement method and the half-power points method are simple, but not 
very accurate. For obtaining more accurate results, either the 
Impulse response method or the Improved resonant forced vibration 
technique should be used. The Improved resonant forced vibration 
technique consists of the estimation of the damping ratio from an 
iterative least-squares curve-fitting procedure for the measured 
frequency response data. The non-resonant forced vibration technique 
Is applicable only for testing very small samples In the low frequency 
range. For studying the variation of damping with frequency, the 
input power-method can be used to get some approximate idea. Finally, 
several new methods developed specially for evaluating the damping of 
viscoelastic materials were described. 
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v- jhe choice of the technique to be used for a particular 

application depends on several factors such as : 1) the modulus of the 
material, 2) the damping of the material, 3) the Instrumentation 
w available, d) the frequency and displacement amplitude range of 

Interest, e) the accuracy required in the estimation of the damping, 
etc. As a rule, while measuring the material damping by any method 
w care must be exercised to minimize the errors due to various factors 

In the data collection, and data analysis process. Some of these 
factors are listed below: 

^ 1) Specimen mounting - mount the specimen so as to avoid external 

losses at the supports, joints etc. A free-free boundary condition is 
ideal for material damping measurements. 

^ 2) Transdu cer/exciter mounting - location and orientation - use of 

mini -accelerometers or non-contacting type of transducers and exciters 
Is recommended to minimize the effects of the added mass on the 

^ response of the system. 

3 ) External distu rbances - Minimize the effects of external 
Influences like, electro-magnetic interference, foundation vibration, 

v acoustic noise, cable motion, air damping, fluid fiow, etc. 

4) No n-linearity of the structure — such as, rattling, banging, 
whipping of loosely fastened components, violation of "small- 

v displacement" theory, non-linear stiffness, non-linear damping, etc. 

I should be avoided wherever possible. 

5) Da ta analysis - The various signal processing errors that could 

v, arise because of inappropriate selection of the measurement parameters 
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such as weighting function, number of averages, amount of overlap, 
calibration factor, etc. should be minimized. Errors due to aliasing 
and leakage should also be reduced. 


III. LITERATURE REVIEW 


Previous Pertin ent Work on the Damping of 
' Composite Material's - 

Efforts to measure, model and Improve the dynamic properties of 
advanced fiber-reinforced polymer composite materials have been 
reported in several publications in recent years. There Is a wide 
range of specimen configurations, materials, frequency range, and 
testing conditions, in the existing literature. Bert [55], Gibson and 
Wilson [56], and Gibson [57] have reviewed the work of various 
investigators in this area. Schultz and Tsai [58] have reported the 
measurement of elastic moduli and damping ratios of glass epoxy 
composite materials. Data on moduli and damping ratios were gathered 
by studying the free and forced transverse vibrations of cantilever 
beam specimens In the frequency range from 5 Hz to 10,000 Hz. The 
elastic moduli showed little tendency to change within this frequency 
range, and the damping ratios were typically of magnitude 10 . They 
however, exhibited both frequency and amplitude dependence. Jones 
[59] has measured the damping ratio of steel epoxy composite material 
having controlled volume fractions and varying wire sizes. 

In a study of dynamic properties of graphite epoxy composites, 
Rehfield et al [60] noticed that these materials experience a 
degradation of mechanical properties due to moisture absorption and 
presence of elevated temperature environment. Leung [61] has studied 



the effects of gamma-radiation exposure on the mechanical properties 
of graphite epoxy composites. Interlaminar shear strength and the 
damping ratio showed an initial increase, followed by a decrease as 
the dosage of the gamma radiation was increased. In another study , 
Springer and his co-workers [62-64] have made an extensive 
investigation on the effect of moisture absorption and desorption on 
the dynamic properties of graphite epoxy and glass epoxy composites. 
The effect of these properties on the damping capacity, however, was 
not included in their study. 

The effect of temperature on the damping and modulus of composite 
materials has also been studied. The variation of Young's modulus 
(E), and loss factor 17 with temperature at fixed frequency and at low 
cy c 1 i c ■ s t ral n amplitude are typically of the form shown In Figure 5 
[65] .- 
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Figure 5. Variation of Storage Modulus and Loss Factor With 
Temperature for Typical Damping Materials 




Three distinct regions are observed, namely the glassy, transition 
and rubbery regions. In the glassy region E is high and rj Is low; in 
the transition region E varies rapidly with temperature and 17 Is high, 
in the rubbery region E varies more slowly with temperature and rj is 
lower than in the transition region, although not always as low as In 
the glassy region. For many materials, not one but several transition 
regions may occur, which usually merge into one another. Yakovlev and 
Bykovskii [ 66 ] have studied the effect of temperature on glass fiber 
epoxy composites. Their results indicate that the variation of 
damping with temperature depends on the kind of fiber, fiber length 
and orientation. The authors also noted that, at cryogenic 
temperatures, the damping of the glass epoxy composites was much lower 

than the damping at room temperature. 

The effect of fiber orientation and lamination on the elastic 
modulus and damping capacity of composites has been studied by some 
Investigators [67-73]. In most composite materials, it has been found 
that the damping is minimum and stiffness is maximum at 0 ° fiber 
orientation, i.e., In the fiber direction. As the orientation angle 
Is Increased, the general trend is for the damping to Increase up to a 
certain value, and the stiffness to decrease. For carbon fiber epoxy 
composites, the maximum damping is reported to occur at a fiber 
orientation of about 67°. According to the study by Mazza et al . 
[ 68 ], cross-ply glass fiber reinforced composites have considerably 
more damping than 0 ° unidirectional fiber-reinforced composites. 
Adams and Bacon [70], in their study noted that the damping of 


graphite fiber-reinforced composites, decreased slightly with the 
cross-ply ratio [ratio of the total thickness of layers of 0° 
orientation to total thickness of layers of 90° orientation]. 

The effect of imperfect Interface bonding on damping 
characteristics of composites has been Investigated by Ki shore, et al . 
[74,75], Very high loss factors were obtained, and the loss factors 
were strongly dependent on the coefficient of friction at the 
Interface. Nelson, et al . [76,77] have also noticed that slip in 
unbonded regions of the interface caused a significant increase in 
damping. 

There is some promise that damping could be used qualitatively to 
detect damage in composite materials since the damping is quite 
sensitive to damage In composites. Gibson and Plunkett [78], in their 
study, progressively damaged cross-ply E-glass composite beams with 
large amplitude vibration and determined changes in damping and 
natural frequency. They found that microstructural damage could 
cause as much as a 350% increase in damping and at the same time 
reduce the natural frequency by 5% or less. Adams [79], has used 
frequency shifts of longitudinal modes in a glass epoxy tube cross 
sectional area. Cawlay [80,81], has extended Adam's work to two 
dimensions, using modal analysis and finite element analysis to detect 
and locate localized damage in cross-ply and single-ply graphite epoxy 
plates. The information on the variation of damping with the number 
of loading cycles has been most inconsistent. The general trend is 
for the damping to increase with the number of loading cycles [82.83], 


In contrast to this general trend, Kim and Matthews [84], noticed a 
rapid decrease In damping In the first few cycles, followed by a 
gradual decrease after that. This Initial rapid decrease In the 
damping was attributed to the residual strain caused by the mismatch 
between the fiber and matrix thermal -expansion coefficients. 

Finally, the mathematical prediction of the dynamic properties of 
laminated composite materials has been Investigated by many authors. 
Hashin [85] has predicted the complex moduli of unidirectional fiber 
reinforced materials based on the properties of resin and fiber. 
Gibson and Plunkett [78] have described a mathematical model for 
predicting the effective complex moduli of unidirectional and cross- 
ply glass epoxy beams in flexural vibration. This model was an 
extension of Hashin' s model. Chang and Bert [86] have Investigated 
the analytical characterization of the damping and stiffness behaviour 
of a single layer of filamentary composite material. 

Adams and Bacon [70], have proposed a model for predicting the 
damping of a narrow angle-ply laminated plate. This method has 
recently been further developed by Ni and Adams [87]. They have shown 
good agreement can be obtained between predicted and measured values 
of damping and stiffness of carbon and glass fiber-reinforced plastic 
beams. In an accompanying paper, Lin, Ni and Adams [88], have used a 
finite element model of Cawley and Adams [89] to predict the natural 
modes and damping capacity of carbon and glass fiber reinforced 
laminated plates. Hwang and Gibson [90,91], have demonstrated the 
application of a finite element model based on the strain energy 



approach for predicting damping In discontinuous fiber composites. 
Johnson and Kienholtz [92], have also used the finite element method 
to predict the damping in structures with constrained viscoelastic 

layers. 

Previous Pertinent Work on the Damping of 
Structural Joints 

It has been found that nearly 90% of the total damping which 
occurs in real structures usually arises in the structural joints 
[93]. The energy dissipation in a joint is a complicated process 
involving several mechanisms, the relative significance of which 
depends on the joint conditions. Because of the complexity of the 
problem it appears that very little work is being done in analyzing 
the damping in joints and built-up structures. In an adhesive bonded 
joint, it is the interface between the adhesive and the adherant that 
Is more prone to environmental degradation and which pi ays .a key role 
in the damping capacity of the joint. A major source of damping In a 
joint is the Coulomb friction damping which arises due to frictional 
forces arising from the relative motion of two contacting surfaces. 
This is usually modelled by a constant force which is proportional to 
the normal load between the surfaces and is directed against the 
velocity vector at each instant. But, in actual practice, the amount 
of energy dissipated depends on both the normal and tangential forces 
in a complicated and non-linear fashion [94-99]. Another source of 
damping in a joint is that due to pumping of air, water or other 
fluids through highly constrained passages [100,101]. The damping is 


caused by the change In spacing of the parts. For example, If a rib 
is riveted or spot-welded to a plate, it will not make contact at all 
points. If the combination is then bent, the clearance between the 
two will change, thus 'pumping' the surrounding fluid (air or water) 
through the narrow passage between them. The viscosity of the fluid 
will cause the damping. 

The damping in a joint is generally non-linear. Most attempts 
are made to linearize the equations in some way [102,103]. A very 
useful method is to calculate an equivalent viscous force component to 
obtain the enegy dissipation. Crawley and co-workers [104-105], have 
developed a method called the force-state mapping for Identifying the 
non-linear properties of structural members such as joints. The 
technique Includes the use of very accurate Instrumentation to measure 
th.e force transmission properties of a joint as a function of the 
relative displacement and velocity across the joint; i.e., as a 
function of the full mechanical state of the joint. The force-state 
map of a general linear spring mass damper sysuem would be a plane 
whose slope with respect to displacement would be the linear stiffness 
K and with respect to velocity would be the linear viscous damping. 
Any deviation in a force-state map from a flat plane is an Indication 
of a nonlinearity in the system. The force-state mapping technique has 
been shown to provide a good method for characterizing the dynamics of 
joints and structural elements. 


IV. MEASUREMENT OF DAMPING OF GRAPHITE EPOXY MATERIAL 


The main objective of the work reported in this chapter was to 
determine experimentally the material damping of the graphite epoxy 
material used in the construction of the truss system of the Hubble 
Space Telescope. Other objectives were: a) to study the effect of 
outgasslng (moisture desorption) on the damping of the material, and 
b) to determine the influence of temperature (both high and low) on 
the material damping values. Damping measurements were made on both 
tube and beam specimens. These samples were supplied by the NASA 
George C. Marshall Space Flight Center. The following sections 
describe the test strategy, experimental apparatus, test procedure and 
results and discussions. 

Tube Specimen 

The composite tube supplied by NASA was in the form of a 
cylindrical tube with 6.17 cm outer diameter, 0.16 cm wall thickness, 
and 95.5cm length. The tube was tested to evaluate Its damping ratio 
value under normal atmospheric conditions and in a vacuum chamber to 
simulate the conditions in space. 

Damping Measurements in Normal Atmospheric Conditions 
Both free vibration and forced vibration methods were employed. 

The ambient temperature averaged 27°C during the whole period of 


experimentation. The tube specimen was first tested with two 
different types of edge support conditions. They were: 1) knife edge 
supports to simulate Idealized simply-supported boundary conditions at 
the edges, 2) a three-points support to hold the tube at the edges. 
Then a free-free support condition with no constraints at the edges, 
was utilized which was found to be most suitable and consistent for 

accomplishing our main objective. 

The idea of using knife edge supports at the two ends of the tube 
specimen was to simulate the Idealized simply supported boundary 
condition. This type of boundary condition Is normally used In the 
theoretical analysis for its mathematical simplicity. But to simulate 
such a boundary condition is rather difficult in practice. The 
constraint here is that the support should not allow any axial or 
vertical deflection of the specimen at the ends, and allow any slope 
at the ends due to any bending moment. With this In mind, it was 
decided to use circular knife edge supports to simulate this boundary 
condition. The design was such that one knife edge would support the 
tube from the inside and the other would support it from the outside 
at each end location. Each circular outer knife edge was made of 
stainless steel In two halves to facilitate easy mounting. Each Inner 
knife edge was also made from stainless steel but was only in one 
piece. A housing to hold the knife edges was designed and fabricated. 
This housing allowed the movement of the knife edges in the axial 
direction for fine adjustment. The housing was made from mild steel. 


Figure 6 shows the knife edges and housing which were designed. Two 
such supports were fabricated, one for each end of the tube. 

Some preliminary measurements of the damping value of the tube 
specimen with the knife edge supports indicated a large contribution 
from the edge supports. Hence, in order to reduce the influence of 
the support damping and to measure only the material damping, a new 
type of edge support was designed and fabricated as shown in Figure 7. 
The idea here was to support the tube at each end at three points 
which were 120 degrees apart. 

Free Vibration Experiments 

A schematic diagram of the experimental set-up used in the free 
'vibration tests is shown in Figure 8. The test set-up consists of a 
non-contacting type (capacitive) displacement transducer (ASP-20) 
along with a measuring amplifier (Accumeasure System 1000 made by MT 
Instruments). A holder was developed for the pick-up to allow for 
easy adjustment of the gap between the specimen and probe surface. 
The output in the form of voltage from the transducer was fed to the 
dual channel signal analyzer (B&K type 2032). The specimen was 
excited for free vibration with an impulse using a small hammer and 
the vibration displacement as a function of time was recorded on the 
analyzer. The damping ratio value of the tube specimen was obtained 
from the time domain data by the logarithmic decrement method. The 
accuracy of estimation of the damping value was later Improved by 
fitting a curve to the data obtained using a least-squares error 
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The experiments were conducted using both the knife edge and 
three-points supports. The specimen was held in the supports at both 
the ends. The tube specimen with the supports was clamped to a rigid 
stand with a heavy base specially fabricated for this purpose. With 
the knife edge supports, three sets of results were obtained, namely 
1) for very large amplitudes of displacement, 2) for moderate 
amplitudes of displacement and 3) for very small amplitudes of 
displacement. The values of damping ratio obtained using the 
logarithmic decrement method In the three cases are summarized in 
Table 1. It can be seen that the damping ratio Is of the order of 
1.0% in all the cases. In order to obtain a better estimation of the 
damping ratio value, a least squares curve-fit was used for the time 
domain data as mentioned earlier. Figures 9 through 11 show the 
least-squares fit for the experimental data obtained with the knife 
edge supports. Note that in Figures 9 through 11 the displacement 
amplitude is plotted on a logarithmic scale. The tube damping ratio 
value obtained from the least-squares fit procedure agrees with that 
obtained from the logarithmic decrement method. It Is seen that in 
Figure 11, for very small displacement amplitudes, there seem to be 
some oscillations in the decay data. This is perhaps due to the 
external noise, since at very small displacement amplitudes, the 



Computation of Damping Ratio of the Tube Specimen From 
Different Cycle Ratios 


Casel 

Case 2 

Cycle 3 

Max. Displacement 

Max. Displacement 

Max. Displacement 

_3 

3.6 x 10 cm. 

1.3 x 10 cm. 

_3 

0.1 x 10 cm. 

r =1.09% 

0-35 

r =1.07% 

0-35 

=1.19% 

0-20 

c =1.00% 

0-20 

C =1.13% 

0-20 

£ =1.47% 

0-10 

r =0.97% 

0-10 

£ =1.16% 
0-10 

=0.91% 

10-20 

£ =1.04% 

10-20 

£ =1.10% 
10-20 

£ =1.09% 

0-15 

=1.20% 

- . 20-35 

£ =0.99% 

20-35 

- 

£ =1.06% 

average 

£ =1.09% 

average 

£ ' =1.16% 
" average 


Legend: the subscript 0-35 indicates that the damping ratio was 

obtained from the vibration decay measurement from 0 to 35 cycles. 
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Figure 9. Free Vibration Decay Curve for the Tube Speci 
(Maximum Displacement = 3.6 x 10~*cm) 
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Figure 10. Free Vibration. Decay Curve for the Tube Specimen 
(Maximum Displacement = 1.3 x 10 ~ J cm) 
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Experiments were repeated on the tube specimen using the three- 
points supports at the ends of the tube. The same experimental set up 
as described earlier was used in this experiment. The damping ratio 
value with this type of support was also about 1.0%. Experiments with 
this type of support showed that the damping ratio value Is very 
sensitive to the positioning of the three supports. For the case of 
the knife edge supports, a Fourier transform of the time domain data 
indicated a dominant resonance peak at a frequency of 390 Hz, when the 
length between the supports was 91.4 cm (36 Inch). For the case of 
the three-points support, the resonance frequency was reduced to 299 
Hz for the same length between the supports. 

The theoretical resonance frequencies of the tube specimen with 
simply-supported ends for flexural vibration were determined using 
beam theory. From this analysis, the first and second modal 
frequencies were found to be 321 Hz and 1287 Hz respectively. Hence, 
it is seen that the experimental resonance frequency of 390 Hz, is not 
matched very closely by that of the theoretical first bending mode 
frequency of .the tube specimen. The discrepancy between the two 
values can be attributed to several factors, such as the knife edges 
not exactly simulating an ideal simply-supported boundary condition, 
the material modulus being different in different directions of the 
tube specimen. 



Forced Vibration Experiments 


A schematic diagram of the experimental set-up used "for 
conducting forced vibration tests is shown In Figure 12. In this set- 
up the specimen was excited by means of a shaker (B & K 4809) driven 
by a random noise generator. The shaker was suspended using a very 
soft rubber rope thus enabling the natural frequency of the shaker 
suspension system to be much less than the natural frequency of the 
test specimen. In the present case the natural frequency of the 
suspension system was observed to be about 1 Hz. The experiments were 
performed on the specimen using the same edge supports as in free 
vibration tests, namely, the knife edge supports, and the three points 
supports. The - frequency response of the specimen in the form of its 
receptance (displacement/force) was measured using an impedance head 
(B&K 8001) and a dual channel signal analyzer (B&K 2032. ). From the 
response curve the half power points and hence their bandwidth of 
frequency separation (Af) and the natural frequency (f;) of the 

specimen were measured. The damping ratio ($•) Is given by £ = Af/2f n . 

The improved half-power points method using a curve-fit analysis 
(chapter 2) was later used in all the cases to obtain a better 
estimation of the damping ratio. The results were not very 
consistent, the lowest value of the damping ratio measured with this 
support was of the order of 0.57%. Using the same experimental set- 
up, a second set of experiments was conducted on the tube specimen 
with the ends supported by means of three-points support. The damping 









ratio value obtained In this case was of the order of 1.15%. The 
difference In the damping ratio of the composite specimen with the two 
edge support conditions confirms that the edge supports had 
considerable influence on the damping ratio of the specimen. 

Having established that the damping ratio of the composite tube 
is strongly dependent on the end support conditions, it was decided to 
support the specimen so that the two ends of the tube were not 
constrained at all. This type of support simulated free-free boundary 
conditions at the two edges. A schematic diagram of the experimental 
set-up for this case is shown in Figure 13. The specimen was 
carefully mounted directly on the shaker and the impedance head using 
a supporting ring (made in two halves) with an Interior knife edge as 
shown in Figure 13. The inner knife edge of the ring allowed the tube 
to be held firmly in the ring. The ring also had a tapped hole on its 
outer surface to connect it to the Impedance head and shaker assembly. 



Figure 13. Tube Specimen with Free-Free Boundary Conditions 




The tube with the free-free edge conditions was then excited 
using a random excitation signal in the frequency range of 0-800 Hz. 
A dominant peak at 509 Hz was observed in the frequency response plot 
as shown in Figure 14. The response was later measured by zooming 
into this peak to get better frequency resolution. The so-called zoom 
analysis is actually a band selectable Fourier analysis, In which 
Fourier transform based digital spectral analysis is performed over a 
frequency band whose upper and lower frequencies are Independently 
selectable. In the baseband analysis, Fourier transform is computed in 
a frequency range from zero frequency to some maximum frequency. This 
digital Fourier transform is spread over a fixed number of frequency 
lines (typically 1024) which therefore limits the frequency resolution 
between lines. Zoom analysis can provide an improvement in frequency 
resolution of more than a factor of 100, as well as a 10 dB Increase 
in dynamic range compared to baseband Fourier analys1s[19] . In the 
present case, the frequency difference between any two consecutive 
measurement points on the analyzer was decreased from 1 Hz to 0.125 
Hz, through the zoom analysis. With this set-up the damping ratio 
value of the specimen was measured and found to be about 0.13%. This 
was the mean value from 10 trials. There was a variation of about 
±0.03% in the values of the damping ratio £. This value is about 8 
times lower than that obtained earlier, and hence was believed to be 
mostly due to the material damping of the specimen as there were no 
structural joints and the specimen was not constrained at its ends. 









The measurements were repeated several times on different days, in 
order to check the consistency of the result. Similar results were 
obtained on all the occasions. The experiment was also repeated with 
two other types of excitation signals, 1) pseudo random noise and 2) 
variable sine with manual sweep. Similar results were obtained with 
both of these types of excitation. 

The next course of investigation was directed at calculating the 
damping induced by the supporting ring used in the above free-free 
experiments. In order to check the damping induced by the supporting 
ring, measurements were made with the top half of the ring removed. 
The specimen then was supported only by the bottom half of the ring. 
It was observed that the natural frequency of the specimen (first 
mode) increased to a value of 552.5 Hz due to the influence of the 
reduced mass of the supporting ring. But the damping ratio value 
obtained was the same as that obtained with the full supporting ring, 
namely 0.13%. This gave more confidence In the measured value of the 
damping ratio of the specimen. The dominant peak at 509 Hz as shown 
in Figure 14 corresponds to the first bending mode of vibration of 
the tube. This was revealed by carrying out a modal analysis of the 
tube specimen. The measured mode shape is shown In Figure 15. 

The next set of experiments was conducted on the cylindrical tube 
specimen using the free-free boundary conditions (half ring) in order 
to study the variation of the damping ratio value with different modes 
of vibration. Figure 16 shows the frequency response curve of the 
specimen in the frequency range of 0 to 3.2 kHz. The numerical value 
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of the resonance frequency of each mode of vibration cannot be 
accurately estimated from this plot due to Its large frequency range 
and hence poor frequency resolution. To reduce this problem, two 
techniques may be used, a) reduce the frequency range of analysis, b) 
zoom the frequency response plot at the natural frequency of interest. 
The resonance frequencies estimated from the plots obtained through 
zoom analysis for the first six modes of the specimen are given in 
Table 2. It is seen from the table that the damping ratio value 
Increases with the mode number. This Increase could be attributed to 
Increased frictional losses at the center support due to Increased 
curvature of the mode shapes at the support point. 


Table 2 

Damping of the Tube Specimen in Different Modes 
for Free-Free End Conditions 

Peak No. Frequency(Hz) Damping Ratio 


1 

552.5 

0.13% 

2 

1514.0 

0.25% 

3 

1873.0 

0.33% 

4 

2204.5 

0.35% 

5 

2232.0 

0.36% 

6 

2841.8 

• 0.36% 
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Figure 16 . Frequency Response of, the T'ube Specimen Showing Several Resonances 



To summarize, all the above experiments were performed on the 
cylindrical graphite epoxy tube specimen under normal atmospheric 
conditions. The damping ratio value was measured using free and 
forced vibration technigues. All of the damping results from the 
above experiments are summarized in Table 3. The damping values 
obtained using the free-free boundary condition for the edges are 
believed to represent most closely the material damping of the 
specimen. As stated earlier the damping value Is significantly 
Influenced by the edge support conditions. Edges and support fixtures 
can normally be expected to Increase the value of damping. 

Tabel 3 

Damping of the Tube Specimen for the First Bending Mode 
Using Different End Conditions 

End Condition Damping Ratio 'Freguency 

Knife Edge Supports 0.90 - 1.10 % 390 Hz 

Three-points Supports 1.00 - 1.20 % 298 Hz, 


Damping Measurements in a Vacuum 

A vacuum chamber as described in this section was designed and 
built to test the cylindrical graphite epoxy specimen. The following 
factors were considered for the design of the chamber: 1) The chamber 
should have a small volume which should thus require a short time to 
create a vacuum In the space. A small volume In the chamber should 
also make it possible to heat the specimens and degas them more 
easily. 2) It should have a minimum number of connecting junctions In 
order to avoid leakage problems. 3) It should be made so that it is 
convenient to excite the cylindrical graphite epoxy specimens With 
different edge support conditions. 4) It should be made convenient to 
test the cylindrical graphite epoxy specimens of different lengths. 5) 
It should be rigid and should be Isolated from the vibration of 
structural parts. Keeping these factors in mind, two designs were 
considered: 1) A rectangular box type enclosure for the specimen, il) 
A cylindrical chamber with end flanges. In view of the difficulties 
that were envisaged in providing vacuum sealing around all the edges 
and also the difficulties that are likely to be encountered in 
providing feed-throughs, a rectangular type of chamber design was 
discarded. 

It was decided to use a cylindrical chamber design whose details 
are shown schematically in Figure 17. It was also felt that it would 
be easy to procure cylindrical stainless steel tubes which would then 
make it possible to avoid many welding and sealing problems. A four- 
way stainless steel cross (15.24 cm in diameter) was used to support 
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Figure 17. Cylindrical Chamber for Testing the Tube Specimen in a Vacuum 



the chamber as shown in Figure 17. The feed-throughs, ports, etc. 
were provided at the top of the cross. A heating coil was provided 
inside the chamber for heating the graphite epoxy specimen. This 
heating coil was necessary to ensure that the composite material 
outgassed rapidly. Two six-inch stainless steel tubes (of 33cm 
length) with end flanges were connected to the flanges on the two 
sides of the cross as shown in the figure. The bottom flange of the 
cross was connected to the vacuum pump unit. The joints at the 
flanges were provided with 0-rings and were sealed with vacuum grease. 
Considerable time and effort was spent to detect the leaks in the 
vacuum system and in the vacuum chamber. The system was able to 

produce a vacuum pressure of about 10 torr (1 torr s 1 mm of Hg) . 

The tube specimen was carefully placed and positioned in the 
vacuum chamber. A free-free edge support condition was simulated for 
the specimen by the use of the full ring with an Inner knife edge as 
discussed earlier. The tube specimen was excited by means of an 
electrodynamic shaker. Figure 18 shows the fixture designed for 
exciting the tube specimen In the chamber for forced vibration 
experiments. The excitation force from the shaker was transmitted to 
the composite tube with the help of a thin rod and a metallic bellows 
assembly. The end flanges of the bellows were connected to the mating 
flanges of the connecting elements (half nipples) by using copper 
gaskets and screw nuts. Care was taken to see that the specimen was 
in a horizontal position and that it did not touch the sides of the 



Figure 18. Experimental Set-up for Exciting the Tube Specimen 
in a Vacuum Chamber 






vacuum chamber. Figure 19 shows a photograph of the vacuum chamber In 
its assembled form, and Figure 20 is a cut-out view of the chamber 
showing the composite tube, heating coil and the Impedance head. 

Experimental Procedure and Results 

As before, the frequency response measured In the present case Is 
also the receptance, which is the ratio of the displacement to force. 
The damping ratio and the undamped resonance frequency were estimated 
using the half-power points method along with the curve-fitting 
technique which was discussed before. Random excitation with a 
Hanning window was used in all the experiments. Some experiments were 
repeated with a swept sine excitation and similar results for the 
damping ratio value were obtained as with random excitation. The 
random excitation was chosen since results are obtained faster with 
th'is type of excitation than with sine sweep testing. The response 
curve was zoomed in the frequency range of 450 Hz to 550 Hz to obtain 
the first bending mode of vibration. This was found to be the best 
set-up for a good frequency resolution. 

First, measurement of the damping ratio of the tube specimen was 
made while the tube was in the vacuum chamber before the air was 
pumped out. Several measurements were made in order to check the 
repeatability of the results. Figure 21 shows the result of a 

measurement made on the tube at 27°C (80°F) under atmospheric 
conditions before the chamber was evacuated. The damping ratio 









measured was about 0.10% and the resonance frequency was about 531.9 
Hz. This value of the damping ratio is somewhat close to the value of 
0.13% obtained earlier for the same specimen in tests performed in 
normal atmospheric conditions. The slight discrepancy between the two 
values can be attributed to the difference In the experimental 
configuration. 

Later the end flanges of the vacuum chamber were closed and the 
vacuum pump was started. The damping ratio value and the resonance 

frequency of the tube specimen were measured on each day for the next 

_2 

two days. The vacuum pressure was observed to be 10 torr on the 
second day and 10" 5 torr on the third day. The temperature inside the 

chamber was still the ambient temperature, namely 27°C (80°F) . The 
damping ratio value of the specimen was measured to be 0.069% on both 
of the days. Figure 22 shows the frequency response plot measured on 
the second day. 

The heating of the specimen was started on the third day. A 
heating coil was placed inside the vacuum chamber below the tube 
specimen (before pumping was started). Care was taken to see that the 
coil did not touch the tube anywhere. The current through the coil 
was gradually increased until the temperature of the space inside the 

chamber was 94^C(200^F). The tube was kept inside the chamber at this 
pressure for a period of one month. The objective was to desorb the 
moisture from the tube material. 
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Figure 22. Frequency Response of the Tube Specimen After Being 
in a Vacuum for Two Days at 27°C 






At the end of the fourth week the heating of the tube specimen 
was stopped. The tube was then gradually cooled to room temperature. 
The damping ratio value and the resonance frequency of the tube 
specimen were measured. When the temperature of the chamber was 

cooled to 26°C (atmospheric conditions) it was interesting to observe 
that the damping ratio value of the tube was 0.063%. This value Is 
very close to the damping value obtained for the tube specimen four 
weeks earlier before the heating was started. 

Finally after cooling, the pressure in the chamber was restored 
to atmospheric conditions and the damping ratio value of the tube was 
measured again. It was found to be 0.083% and the resonance frequency 
was 531.20 Hz. Earlier, under atmospheric pressure, the tube damping 

ratio value was measured to be 0.10%. 

The tube was then carefully dismantled and it was weighed in a 
precision balance. The difference in weight of the tube before and 
after heating was found to be 4.20 grams. The percentage change in 
the weight of the tube from its original weight was 0.9%. This was 
assumed to be caused by the loss of moisture which occurred while the 
tube was degassing. The change in damping ratio value due to 
outgassing of the tube specimen was observed to be 0.016%. 

All of the results discussed so far in this section are 


summarized in Table 4. 


Summary of the Damping Data of the Tube Specimen In a Vacuum 


activity 

CHAMBER 

PRESSURE 

DAMPING 
RATIO • 

MASS OF 
THE TUBE 

1. 8efore pumping the air out 

Atmospheric 

0.10 X 

466.8 grama 

2. Immediately after pumping 

-2 

1 0 lorr 

0.07 X 


3. Continued pumping-2 days 

-j 

1 0 torn 

0.07 % 

— 

4. Continued pumping and 
heating the tube to about 
WO C for 30 days 

-2 -J 
W - W torr 

0,06—0.09 Z 

- 

5. Stopped heating and 
pumping. Atmospheric 
conditions restored 

Atmospheric 

0.08 X 

462.6 grams 


• Mean value of 1*0 trials 1 torr™ t mm Hg 




Temperature Effects on the Damping of the 
Graphite Epoxy Material 

The next course of the investigation was directed at finding the 
effects of temperature on the damping of the graphite epoxy material. 
Experiments were conducted on both tube specimen and some beam-like 
samples using specially built temperature chambers. Forced vibration 
experiments were conducted on these specimens with free-free boundary 
conditions. All of the measurements were made in air, since, It was 
more convenient to conduct experiments at controlled temperatures in 
air than in a vacuum. 


High Temperature Experiments 

Figure 23 shows a schematic diagram of the experimental apparatus 
used to study the effects of high temperatures on the damping of beam 
samples. The specimens were carefully mounted inside the temperature 
chamber and were supported at the center by a steel rod. The steel 
rod was tapped and attached through a hole in the specimen by means of 
a small nut. The other end of the steel rod was connected to the 
shaker via an impedance head. The force and the acceleration signals 
from the impedance head were fed to the FFT analyzer after 
amplification by the charge amplifiers. The temperature Inside the 
chamber was monitored and controlled precisely by two chromel-alumel 
thermocouples and a temperature control programmer. As before, damping 
ratio values were extracted from the receptance plots using the 


improved half-power points method. 








Figures 24 and 25 show the variation of the damping ratio with 
temperature for graphite epoxy beam specimens of two different 
lengths. The temperature where the first peak appears is believed to 
be the glass transition temperature of the material. As can be seen 
from these figures, the variation of £ with temperature has almost the 
same trend in both cases, and also it agrees well with the standard 
plot found in the literature [see Figure 5 ]. 

A similar experiment was conducted on the tube specimen. Because 
of the limitations of the heating equipment, It was possible to raise 

the temperature of the tube only up to 300°C. For better temperature 
control , the temperature of the tube was measured at three points 
along its length. The temperature of the tube was gradually Increased 
from atmospheric temperature and damping measurements were made at 
regular Intervals. The frequency response was measured when the tube 
temperature was the same at all the three points. Figure 26 shows 
the final result for the variation of £ from atmospheric temperature 

to about 300°C. This is again for the first bending mode of vibration 
of the tube. It is seen from this figure that the trend in the 
variation of £ is somewhat similar to that of the beam results, i.e. 
the damping ratio increases with temperature. The discrepancy In the 
absolute value of the damping ratios and the value of the glass 
transition temperature as observed in the three plots could be 
attributed to a) difference in the geometry of the specimens, b) 
difference in the frequencies of vibration and c) difference in the 
lamination (fiber winding) of the beam and the tube samples. 
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In order to gain confidence In the experimental methodology 
followed, and to validate the results obtained so far, It was decided 
to conduct the same temperature experiments on some metal beams. 
Figures 27 and 28 show the variation of damping with temperature for 
an aluminum beam and a stainless steel beam respectively. As can be 
seen from these results, the damping capacity of these materials is 
almost constant with temperature. 

Low Temperature Experiments 

The experimental apparatus used previously for conducting high 
temperature experiments was not suitable for conducting low 
temperature experiments. Hence, an entirely new temperature chamber 
was designed and fabricated. It was initially planned to conduct the 
low temperature experiments in a vacuum chamber. Thus, the set-up 
shown in Figure 29 was developed. As can be seen from the figure, the 
specimen to be tested was enclosed in a temperature chamber, which, in 
turn, was placed in a vacuum chamber. The specimen was connected to 
the shaker (the shaker was kept outside the temperature chamber, but 
inside the vacuum chamber) by means of a steel rod and an impedance 
head as was done during the high temperature experiments. The 
temperature chamber consisted of double-layered side walls forming a 
closed cavity. The cavity had an inlet and an outlet for the flow of 
liquid nitrogen through it. The outside surface of the chamber was 
well-insulated so that only the specimen which is inside the chamber 
will be cooled. There was also a small heating coil in the chamber to 
heat the specimen if required. The flow of liquid nitrogen was 







controlled by varying the pressure In the liquid nitrogen dewar. Dry 
nitrogen gas from a cylinder was sent into the dewar, Increasing the 
pressure and thus forcing the liquid to flow out through the siphon 
Into the chamber. The temperature In the chamber was controlled by 
adjusting the flow rate and also when required using the heating coil. 
All the hoses and tubes were insulated to minimize thermal losses. 
Efforts to maintain a vacuum and low temperature conditions inside the 
chamber were unsuccessful because of many experimental difficulties. 

Hence all the experiments were conducted under atmospheric 
pressure. Experiments were done on the following four samples: 

1) Graphite epoxy beam of dimensions 16.5 x 2.5 x 0.3 cm, 2) Graphite 
epoxy beam of dimensions 25.4 x 2.5 x 0.3 cm, 3) Aluminum beam of 
dimensions 25.4 x 2.5 x 0.3 cm and 4) Graphite epoxy beam of 
dimensions as in 1, with a simple lap joint having an overlap of 8.9 
cm. 

The procedure used in the experiment was as follows. The 
specimen was first carefully mounted inside the temperature chamber 
and was supported at the center by a steel rod as explained before. 
The force and acceleration signals from the impedance head were fed to 
the FFT Analyzer. First, damping measurements were made at atmospheric 
temperature. Then the sample was cooled to the desired temperature by 
feeding an appropriate amount of liquid nitrogen into the cooling 
chamber from the dewar. The temperature of the sample was monitored 
precisely by means of two chromel-alumel thermocouples attached at two 
different locations on the sample. The difference between the 



temperatures at the ends of the sample never exceeded 1 to 2°C during 
the experiment. The frequency response in terms of the receptance 
function (displacement/force) was measured on the FFT analyzer. 
Random signal excitation was used and 100 averages were taken for each 
measurement. The whole measurement process took only about 50 
seconds. But more time was spent in achieving a stable temperature, 
since the temperature control was done manually by adjusting the rate 
of flow of liquid nitrogen. The experiments were quite consistent and 
the results were repeatable. The following paragraph describes the 
results. 

Figures 30 and 31 show sample frequency response curves for the 

second specimen tested at temperatures of 24°C and -130°C 
respectively. The peak at 262 Hz corresponds to the first bending 
mode of vibration and the peak at 1330 Hz corresponds to the third 
bending mode. All the damping measurements were focused on these 
modes. It was not possible to excite even numbered modes since the 
beam was excited at its center which is a node point for even numbered 
modes. 

Figures 32 through 36 show the variation of damping ratio with 
temperature for the four samples mentioned before. The curve-fitting 
method based on the half-power points method was used in each case to 
compute the damping ratio from the receptance plots. Figure 32 
corresponds to the first graphite epoxy beam specimen. The damping 

ratio value of the sample was 0.23% at 24°C and was reduced to about 



half this value (0.12%) when the temperature was reduced to about - 

60°C. There was a small difference of 3 Hz In the modal frequency 
between the two readings. A somewhat similar trend can be found in 
the other plots as well. It is particularly clear from Figure 36 that 
there is a slow and gradual decrease In the damping ratio value from 
an atmospheric temperature to cryogenic temperatures. This is in 
agreement with some data available in the 1 iterature[66] . 
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Figure 30. Frequency Response of the Graphite Epoxy Beam at 
Atmospheric Temperature 
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Figure 31. Frequency Response of the Graphite Epoxy Beam at -130° C 






Damp 1 ng Ratio 


° c 


D3 

c C 


Temperature - uog ^ 

Figure 32. Plot of Damping Ratio of Graphite Epoxy Beam 1 vs. 
Temperatures Below Ambient 
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Figure 33. Plot of Damping Ratio of Gr *P h l te J poxy Beam 2 VS 
Temperatures Below Ambient (Mode 1) 
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Fiaure 34. Plot 'of Damping Ratio of Graphite Epoxy Beam 2 vs 
Temperatures Below Ambient (Mode 3) 
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Figure 35. Plot of Damping Ratio of Aluminum Beam vs. 
Temperatures Below Ambient 
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V. DAMPING OF BONDED STRUCTURAL JOINTS 


Structural adhesive bonding of composites to composites, as well 
as composites to metallic components, has developed rapidly due to 
advances in composite materials and adhesive bonding techniques. 
Structural composite bonded joints are primarily of the overlap type 
(single or double overlap). Scarf joints are structurally efficient, 
but are difficult or costly to manufacture. Butt joints although 
simple to process, are not considered for primary structural joints. 
This is because of their limited load carrying capacity due to the 
presence of cleavage stresses when the loading is not truly tensile. 
Angle joints are not usually preferred but are sometimes necessary. 
Tube joints are difficult and costly to manufacture. 

The paper by Goland and Relssner [106] Is regarded as a classic 
work in the area of static analysis of a simple lap joint. Goland and 
Reissner studied the stresses in bonded single lap joints for two 
different cases. In the first case, the bond layer was very thin and 
had no contribution to the joint flexibility. On the other hand, in 
the second case, the bond layer was so thick that it was the primary 
contributor to joint flexibility. In both cases, they derived 
equations for evaluating the shearing and normal stresses In the bond 
layer as well as those in the jointed plates. They found for equal 
thickness isotropic plates, that the bond layer shear stress has a 
nearly uniform distribution except for large concentrations near the 



ends of the joint. The peel stress (stress perpendicular to the bond 
layer) also had high values near the joint edges although not nearly 
as high as the Inflexible bond case. 

In the Goland and Reissner analysis, the peel and shear stresses 
were assumed constant across the adhesive thickness. In later works 
by Ojalvo and Eidinoff [107], Carpentor [108], Kline [109], attempts 
are made to Incorporate a linear variation of these stresses across 
the thickness of the adhesive. Delale and Erdogan [110] have carried 
out the stress analysis of a bonded lap joint system assuming that the 
adherents are elastic and the adhesive is linearly viscoelastic. 
Renton and Vinson [111], and Delale, Erdogan, and Aydinoglu [112] have 
attempted to include anisotropic adherents in the mathematical model. 

Hart-Smith [113-116] was the first Investigator to advocate the 
use of continuum mechanics In the analysis of bonded joints. He has 
analyzed double-lap, single-lap, scarf, stepped-lap and tapered-lap 
configurations. Tensile, compressive and In-plane shear stresses in 
the system were considered based on an elastic-plastic analysis of the 
configuration. The author has also discussed joint efficiency and 
potential failure modes for each of the above configurations. 

The above review Indicates that much of the work done in this 
area has been confined to the static analysis of bonded joints. 
However, a lot of work has been done in the area of dynamics of 
sandwich beams and plates, consisting of a viscoelastic core material 
constrained between two elastic layers. The reader is referred to 
excellent review articles by Mead [117] and Nakra [118], for more 


Information on this subject. Only recently, two papers have been 
published on the subject of vibrations of bonded joints. The one by 
Saito and Tanl [119], deals with the vibrations of a bonded lap joint 
system and the other by Prucz [120] deals with a quasi-static analysis 
of a symmetrical double-lap configuration. The analytical model 
described by Prucz is similar to the one-dimensional model of Hart- 
Smith for the analysis of a fully elastic double lap joint. The 
author, however, has incorporated the viscoelastic behavior of the 
adhesive layers in the joint. A quasi-static analysis of constrained 
layer damping treatment Is utilized to evaluate the joint damping 
properties. Saito and Tanl [119] have derived equations for 
predicting the modal parameters of the coupled longitudinal and 
flexural vibrations of a system consisting of a pair of elastic beams 
Tap-jointed over a certain length by a viscoelastic material. 
Numerical results are presented for the case of fixed-fixed boundary 
conditions at the ends. A simplified version of this model to Include 
composite beams has been developed in the present investigation. The 
details of this analysis together with numerical and experimental 
results are presented in this chapter. Some experimental results on 
the damping of other types of bonded joints are also included in this 
chapter. 



Analytical Modeling of the Fl exural Vibrations 
of a Bonded Lap-Joint System 

The system chosen for study is as shown in Figure 37. The 
objective is to arrive at a suitable mathematical model to evaluate 
the damping ratios and resonance frequencies of the system for free 
flexural vibration. The system consists of a pair of rectangular 
beams lap-jointed over a certain length. The bond between the two 
beams is achieved by means of an adhesive whose thickness is small 
compared with those of the beams. The unjointed ends of the beams may 
have any physically realizable boundary conditions, but, in this case 
are assumed to be simply-supported. Equations of motion are, first, 
derived assuming the beams to be made of composite materials using the 
complex modulus approach. The adhesive is assumed to be linearly 
viscoelastic and the widely used Kel vin-Voight model is used to 
represent the viscoelastic behavior of the adhesive. Both the complex 
modulus approach and the Kelvin Voight solid model approach are 
Identical In that the constitutive relationship between stress and 
strain is similar to Hooke's law, but Includes a complex rather than a 
real material constant. According to this approach, the deformation 
field induced in a material by a simple harmonic oscillating load 
contains an elastic or storage component, and a dissipative or loss 
component. This model is restricted to cyclic oscillations, and is 
believed to be closely associated with a simple physical 
interpretation of a viscoelastic behavior. 
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The starting point for the development of this model Is the 
consideration of the dynamic equilibrium equations of the overlap 
(joint) region. The equations of motion in the joint region are 
derived using a differential element approach. The transverse 
displacements of the upper and the lower beams are considered to be 
different. The normal force between each beam and the adhesive layer 
is obtained from the Kel vin-Volght viscoelastic model of the adhesive. 
The shear force at the interface between the adhesive and the beam is 
obtained from the simple bending motion equations of the two beams. 
The resulting equations of motion are combined with the equations of 
transverse vibrations of the beams in the unjointed regions. These 
are later solved as a boundary value problem using a knowledge of the 
boundary conditions at the unjointed ends and motion continuity 
equations at the ends of the overlap. The eigenvalues and the 
eigenvectors of the system are obtained numerically by an iterative 
technique using a computer. 

Derivation of the Mathematical Model 
The system is hypothetically divided into three parts as shown in 
Figure 38. The coordinate system chosen for each part is also shown in 
the same figure. 

Part 1 - Overlap Region 

Although, this portion of the system is identical to a three- 
layer sandwich beam, the analysis here is somewhat different from that 
of a three-layer sandwich beam. This is because, in the vibration 



analysis of a three-layer sandwich beam, the displacements of the 
upper and lower beams are usually assumed to be Identical. In the 
present case, however, the system configuration does not permit the 
Inclusion of such an assumption. Several other assumptions are, 
however, made in the present analysis as follows: a) the analysis Is 
carried out only for the free flexural vibrations of the system; the 
longitudinal and rotary inertia effects are neglected, b) there is no 
slip at the bond interfaces, c) all points on a plane normal to the 
longitudinal axis of the beam have the same transverse displacement, 
d) the longitudinal force In the adhesive layer may be neglected 
relative to those forces In the beams, e) since the thickness of the 
adhesive layer is small relative to the thickness of the beam, the 
mass (inertia) of the adhesive may be neglected, and the normal and 
shear forces In the adhesive are assumed to be constant. 

A free-body diagram of a differential element of length dx of the 
composite three-layer part of the system Is shown in Figure 39. N is 
the axial force in the beam, V represents the shear force in the 
beams, M is the bending moment in the beam, s is the shear force per 
unit length at the bond interfaces, and p is the transverse force per 
unit length between the beam and the adhesive. The above quantities 
with subscript 1 refer to the upper beam and subscript 2 refer to the 
lower beam. Furthermore, yj = y 1 (x,t) Is the transverse displacement 

of the upper beam, and y 2 = y 2 (x,t) is that of the lower beam. 



ANALYSIS 



Figure 38. Hypothetical Division of the System Showing the Coordinates 
For Analysis 


pari 1. FBD 



. Free Body Diagram of a Differential Element in the Overlap Region 


Figure 39 
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The sum of forces In the x-direction on beams 1 and 2 yields, 

5N i / x 

Nj + gY - dx - Nj + sdx = 0, and (5.1) 


6Nj 

N 2 + dx - N 2 - sdx = 0. 


(5.2) 


Equating the sum of forces In the y-direction to the corresponding 
Inertia forces, we get for beam 1, 




a j yi 


Vj + Vj + dx + pdx - /ijAjdx g.j. 2 , 


(5.3) 


and for beam 2 


5V ’ „ 

-Vi + v * + aF dx 


0 J y a 

pdx * p 2 Ajdx gp- 


(5.4) 


Fhe term pk in Equations (5.3) and (5.4) Is the mass per unit length, 
rfith p representing the density and A representing the cross sectional 
area of the beam. Summing the moments about the center of the right 
edge, for beam 1, 


3M X h x 

M x " (Mi + dx ) “ sdx 2 + V i dx = 0 » 


(5.5) 


for beam 2, 

QM- h 2 , x 

M 2 - (M 2 + g^- dx) - sdx 2 + V 2 dx = 0 . (5.6; 

The simplified versions of Equations (5.3) and (5.4) are: 

dV i d *yi ' (5.7) 

9x + P = Pi A i at 2 ' 

av 2 32 y 2 

ax - p = Pi*! at*" • 


(5.8) 


Differentiating Equations (5.5) and (5.6) with respect to x, and 
simplifying we obtain, 


avj a 2 Mj \ as 

ax = ax 5 * + T ax ' 

and 

av 2 a 2 Rj h 2 as 

ex = ax 5 " + 2 ax ■ 

From Equation (5.7) through (5.10), we have, 


(5.9) 


(5.10) 


a 2 M 1 hi as a a y x 

ax 2 " + 2 ex + p = PlAl at 2 " * ^ 5 ‘ n ^ 

and 

a 2 ^ h 2 as a 2 y 2 

ax 2 " + 2 ax " p = ^ aAj at 2- • ( 5 . 12) 

Equations (5.11) and (5.12) represent the equations of motion of the 
upper and lower beams respectively. From the classical theory of pure 
bending of beams, the bending moments and Nt, may be related to the 

corresponding transverse displacements, y 1 and y a through the flexural 

rigidity term El, where E is the Young's modulus of the material of 
the beam, and I is the second moment of the cross-sectional area with 
respect to the centroidal axis. Noting that the y-axls is positive 
downwards, from beam theory, 

3 2 yi 3 2 y a 

Mj - - Ej I x and M a = I a . (5.13) 


It should be noted that Equation (5.13) is valid for elastic, 
homogeneous beams having a constant area of cross-section. It turns 



out, however, that an equation similar to the above can be written for 
the composite beams of the present problem, by simply replacing the 
real modulus by its complex counterpart. This analysis is valid only 
for harmonic time dependence [10]. With this in mind, we can write, 

Mj = - E* Ij ^ , and M, = - E^ I, ^ 2 , (5.14) 

dx 2 dx 2 

where Ej = Ej (1+iaw), and Bj = E^l+i/fo/). (5.15) 

a and /? are some constants to be assumed depending on the material, 

and ui is the frequency, and 1 = \M. The next objective is to find 
suitable expressions for p and s In terms of y 2 and y 2 , so that 

Equations (5.11) and (5.12) can be solved. 

Expression for p 

As mentioned before, the adhesive layer Is modeled as' a Kelvin- 
Volght viscoelastic solid, because of the mathematical simplicity. A 
detailed derivation of a standard-linear model for a viscoelastic 
material and its reduction to a Kel vin-Voight model is shown in 
Appendix B. Accordingly, we can write, 

dy 1 9y 2 

P = K (yj -y 2 ) + c (^ - at ) • ‘ ( 5 ' 15 ) 

K refers to the transverse stiffness per unit length, and c is the 
viscous damping coefficient. K may be related to the storage modulus 
E c , the width b, and the thickness h c of the adhesive by, K = E c b/h c . 
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Expression for s 

The interfacial shear force s may be found by considering a small 
element of the viscoelastic layer as shown in Figure 40. 

The element is deformed due to the axial displacements, Uj and u 2 


of the bonded surfaces of beams 1 and 2 respectively, caused by their 
bending moments. For small displacements, the shear stress r, is 
related to shear strain in the adhesive by: 


r = G c { 


“i- U 2 


a (uj- u a ) 
} + ? at "FT 


(5.17) 


The expression on the right hand side excluding G c Is the shear 

strain. 6 is the real part of the shear modulus of the adhesive 

material. As before, 7 is another constant which is included here to 
take care of the viscoelastic nature of the adhesive. 

Assume that the axial deformations, u 2 and u 2 are caused only by 

• the bending of the beams and that there are no external axial forces 

In the beams. Then, from the classical theory of pure bending of 

beams, it can be shown that, 

hi 8 y x , h J a >2 (5.18) 

Ul= 2 dx ' and Uj “ “ 2 9x 

7 

I 
I 

/ 

I 

„_/ 

Figure 40. A Small Element of the Viscoelstic Layer 


Mjl 
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In deriving the above expressions, It is also assumed that the neutral 
axis of the upper beam always remains above that of the lower beam. 




Then, 


hj 3y x h 2 3y 2 
u i‘ u 2 = 2 9x + 2 9x 


(5.19) 


From Equations (5.17) and (5.19), the interfacial shear force per unit 
length s is given by 


G b 


0yi 




3yi 


3y a 


s = 2F [ ( hi ax + h 2 3x ^ + 7 at ^ hl 8x + ha 3x ( 5 * 20 ^ 

c 


It Is interesting to note here that for a simple harmonic time 
dependence of the form, e 

(5.20) can be written in a more familiar form as 


for y 2 and y 2 , Equations (5.16) and 



P = K* [y x -y 2 ] , and 

(5.21) 


G*b dy a dy 2 



s = 2F [h i d7 
c 

(5.22) 


' where, K c = K(1 + i ^ ), and G c = G c (1 + i'ju) . 

(5.23) 


K c and G c may be considered as the complex stiffness and complex shear 

modulus of the adhesive material respectively. 

The next step is to obtain the final form of equations of motion 
(5.11) and (5.12) by making use of the relations developed so far for 

Mj , M 2 , p and s. This is easily done by noting that y x = Yj e^ wt , and 
y 2 = y 2 e 1wt , where Y t and Y 2 are now functions of x only, w is the 


W 
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complex natural frequency, the real part w R representing the damped 
natural frequency and the ratio of the Imaginary part Uj to 


representing the modal loss factor rj. The modal damping ratio £ Is 
simply equal to rj/2. Now, from Equations (5.11), (5.14), (5.21) and 
(5.22), we obtain for the upper beam: 


* 

E i dx-* 


G*bh. d 2 Yj 

[h x - TTT + h 


d 2 Y, 


4h 


dx 2 


2 dx 2 


’] + K* (Yx-Y,) - p, M’Yx = 0, 


(5.24) 


and, similarly from Equations (5.12), (5.14), (5.21) and (5.22), we 
get for the lower beam: 


^ * 2 dx 4 


G c bh > 

4h„ 


d 2 Y x 


d 2 Y, 


Chx dx 2- + hj Hx 2 " 3- K c ^- Y2) - = °' 


- (5.25) 


Here, Y x = Y x (x) , Y 2 = Y 2 (x) . 

Equations (5.24) and (5.25) are coupled equations of motion of the 

system, the solution of which can be obtained by assuming a solution 

Lx X n x 

of the form Y x = A n e , and Y 2 = B n e 


Non-dimensional jzati on 

The above equations of motion (5.24) and (5.25), should be 
expressed in a non-dimensional form so as to avoid overflow problems 
on the computer during the solution scheme. This is done in this 


section for the special case in which the two beams are assumed to be 
Identical . 

Let E x = E 2 = E, and h x = h 2 = h. The following non-dimensional 
parameters are used: 


- Y - Y - x _ 0 

Yrr'l V = 2 Y= 0 = ^ 

V ‘ 5 ' 1 E L ' L ' 


(5.26) 


where L = H + H = length of the beam and Z c = length of the 

adhesive layer. 

Also, let 


- G c - E c - h - ^c 
G c"E ,E c = E ’ h " L ’ £ c = L * 

The non-dimensional natural frequency is expressed as 


(5.27) 


El 


u = - , where w 2 = 


(5.28) 


-o 0 Mo- 

using, Equations (5.26) through (5.28), Equations (5.24) and (5.25) 
written in non-dimensional form are: 


d<Yx 

dx 4 


_ * 
3G c 


h F c dx 2 


d 2 Y. d 2 Y 2 12 E_ 

[ -i- + — ] + — £ (Yj 


- 3 + = r 

dx 2 h c K 3 


- Y 2 ) - w 2 Y I =0, (5.29) 


and 


d 4 Y, 

3G* d 2 Y. d 2 Yj 

_* 
12 E r 


1 

=r£ [ — + — 1 

(Y, - Y- ) - w 2 Y 2 = 0, 

(5.30) 

dx 4 

h H c dx 2 dx 2 

=T| 

O 

i 


_* 

E c (l+iJ? a ) _* 

G c ( 1+1 fc) 


Here, E c 

= eTTTT^TT ' ind G c 

" ECl+1^1 ) ’ 

(5.31) 
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^ may be considered as the modal loss factor of the beam material, 
and 7 j 2 and tj 3 as the modal loss factors of the adhesive material in 
bending and shearing motions respectively. rj l , ifo and rj t may be 

assumed suitably depending on the material, while computing the 
eigenvalues and eigenvectors for each mode. 

As mentioned before, Equations (5.29) and (5.30) can be decoupled 
by a series solution of the form 


Ax _ Xx 

Y x = A n e n , and Y a = B n e 


(5.32) 


where A and B are constants to be determined from the boundary 
n n 

conditions. 

Substituting Equation (5.32) in (5.29) and (5.30), we obtain, 


3G* 12E* 

Xi A_ — [ X* A„ + XJ, B n ] - — ! = [B„-A n ] - * A„ - 0 (5.33) 


n n 


Tf H, 


h_ h 3 

c 


* — * 

3G 12E 

X< X* A n ♦ X«„ B n ] - — £ [B„-A n ] - V B„ - o 

K F c he h» 

The above equations written in a matrix form: 


(5.34) 


(X< " C i X n + C s ‘ w2 ) ~( c i +c s) 

-(Cj+ C f ) (XJ - C.XJ + c 2 - *) 




(5.35) 
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/< 
i ♦ 


V' 


where, C t 




(5.36) 


For a non-trivial solution, the determinant of the above matrix set 
equal to zero yields eight roots of X, i.e., 

(X< - CjX* + C 2 - w*) S - (C a + C 2 ) J = 0. (5.37) 

This can be simplified as 

(X« " C i X n + 2 Cj + C i " w 2 ) =0 (5.38) 

and 

(A< ' C i X n ~ ~ c i) 3 0 . (5.39) 

Equations (5.38) and (5.39) each yield four unique roots for X. 

Furthermore, from the above matrix, the constants A„ and B„ are 

n n 

related by: 

B n = h A n. 

n r n n f 

where, 


6 = [ 
r n L 


K - C . X ’n + C > 


c + c 

'-l ^2 



Finally, we have 


Yi 


8 

-E 

n=l 


V 


(5.40) 


(5.41) 


and 
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V 


v, - E f>A e V . (5.42) 
n=l 

Equations (5.41) and (5.42) now have only eight constants to be 
determined by applying the boundary and continuity conditions as 
described later in this chapter. 


Analysis of Parts 2 and 3 

The differential equations of motion for the transverse vibration 
of the beam portions of the system as shown in Figure 38 are: 


9<y 


J 8 *j 


M + p A — = 

^j A j at* 


j = 1,2. 


(5.43) 


For the special case of identical beams, 

Dj = D 2 = E*I , and p x A 1 = /> 2 A, = pk. 

Assuming y bj = Y bj (x) e 1wt , j=l,2, Equation (5.43) becomes 

d<Y h , pk 

asr* - W Y bj - J ^ 

vJ 


(5.44) 


In non-dimensional form, 


^ - u 2 Y M = 0, j = 1,2 . (5.45) 

dxij DJ 

J 

The solutions of the above equations are obtained, as before by 
assuming a series solution of the form: 
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Y bl 'Jj C n ^ • and Y b2 D n ^ ■ (5 ' 46) 

where 6 * = w 2 , and C and D are constants to be found from the 
“n n n 

boundary conditions. 

Equations (5.41), (5.42), and (5.46) have 16 unknown constants. The 
following boundary and continuity conditions are applied to determine 
those constants and to obtain the frequency equation. 

Boundary Conditions 

Considering simply-supported boundaries at the un-jolnted ends of 
the two beams, at points 1 and 2 (see Figure 41 ), the transverse 
displacement and bending moment are zero. 


*bi 

> 

\ 

7 


x„ ' 

- 


J 

> 

fY 2 > 

w n 


Figure 41. Illustration of the Locations of Boundary and Continuity 
Condi tions 
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d 2 Y.. 

1 .e. , a) at x 1 = 0, Y.. = 0, and — = 0, (5.47) 

U1 dxf 

d 2 Y h? 

b) at x 2 = l, Y b2 = 0, and — 0. (5.48) 

d Xj 

Next, points 3 and 4 are free ends. Hence, the bending moments there 
are zero. Furthermore, the shear force at these points can be 
obtained from Equation (5.22). Hence, we obtain, 

c) at x = 0 (point 3) , 


d 2 Y, 



d 3 Y a 

_* 

6G r 

dY x 

dYj 

(5.49) 

2 

— 

0, and 


c 


+ — 1 , 

dx 2 



dx 3 

" H 2 

dx 

dx 


(point 

4), 






d 2 Y. 



d 3 Y, 

6G r 

dYj 

dY. 


i 

= 

0, and 

1 

c 

[ -zr 

+ — ] • 

(5.50) 

dx 2 



dx 3 

" F 2 

dx 

dx 



Continuity Conditions 

Referring to Figure 41, at points A and B, we equate the 
y displacements, slopes, bending moments, and shear forces for the two 

/; hypothetical sections of the beams. These are mathematically 

• I 

expressed as: 


= a and x = 0 (point A), 




dY1 

. dY bi 

dx 

dXj 

d 2 Yj 

. d,Y bi 

dx 2 

dx; 

d 3 Yj 

. d ’ Y bi 

dx 3 

dxjf 

= 0 

and x » 

II 

Y b2 ' 

d 

. dY b2 

dx 

dx 2 


(point B) 


( 5 . 51 ) 


( 5 . 52 ) 
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Equations (5.47) through (5.52) yield 16 equations containing the 16 
unknowns, viz., A lf A,, . . . A,, C lt C s . . . C 4 , and D 1( D, t . . D 4 . 

These equations can be expressed In matrix form as: 



The elements d.. of the matrix D can be written in groups as follows: 
• J 

1) d.. j » 0 for 1 = 1 ... 4 and j = 1 ... 8 

i = 1,2 and j = 13 ... 16 

i = 3,4 and j =9 ... 12 
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1 =5 ... 8 and j =9 ... 16 
1 =9 ... 12 and j = 13 ... 16 
1 = 13 ... 16 and j = 9 ... 12 
2) d. . = 1 for 1=1 and j = 9 ... 12 

* J 

for 1 = 9 and j = 1 ... 8 
d tJ = -1 for 1 = 13 and j = 13 ... 16 


3) d 2J = /3J_ 8 for j = 9 ... 12 



P* 12 2 - 

4) d 3j = e J for j = 13 ... 16 

5) d 4j = -12 d 3j for J = 13 ... 16 

6) d 5j . = for j = 1 ... 8 

„* 

6G 

7) d 6j ' ^ < 1+ W for j = 1 ... 8 
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8 ) d 7j = X J' 6 J f0r J = 1 *•* 8 

6G* U c 

9) d gJ . = [XJ - — (1 + f,)*j]e for j 

10) d 10j = x j ^ or J = 1 ••• 8 

11) d^j » XJ for j = 1 ... 8 

12) d 12j . = XJ for j = 1 ... 8 

X 1^c 

13) d 13j . = J for j = 1 ... 8 
H) d HJ = Xjd 13J for J = 1 ... 8 

15) d 15J - XJ d nJ for j • 1 ... 8 

XX 

16) d i6j ’ e J d fj for J » 1 ... S 



125 


17) d 9j ’ ' e ^ J " 8 ^ for J " 9 ••• 12 

18) d 10J - Aj_gdgj for J - 9 ... 12 

19) d nJ - ^.. 8 d 9J for j ■ 9 ... 12 

20) = ^j-8 d 9j ^ 0r ^ ^ *** ^ 

2!) d uj = - Pj_ 12 for j = 13 ... 16 

22) dj^j = - ^ or j = 13 ... 16 

23) d 16j = ' ^j-12 for J = 13 ••• 16 

The determinant of the above matrix D, set equal to zero gives the 
frequency equation, 1 . e . , for a non-trlvial solution, 

0et [D] 16x16 * * °‘ (5 ‘ 54) 

The roots of Equation (5.54) yield the complex natural frequencies 
(w n s) of the system. 

*nj = “RJ + 1 “Ij- 

where j = 1, 2, . . . represents the mode number. 



= modal resonance frequency, and 


w Rj 

W Ij 

n. = = modal loss factor. 

w Rj 

Once the eigenvalues are evaluated, the corresponding eigenvectors 
needed to plot the mode shapes can be found by using any 15 of 
Equations (5.53) with one arbitrary constant. 

Solution Procedure 

The solution of Equation (5.54) was obtained numerically on a computer 
using an iteration method. A computer program has been developed on a 
Harris 800 system to obtain numerical solutions. 

The complex eigenvalues (w R + 1 u^) are solutions of the complex 
equation ' 

D r (w r + 1w.) + 1 Dj(w r + 1 ^0 = 0, (5.55) 

where D d and D T are the real and imaginary parts of the determinant of 
K 1 

the matrix D. The inputs to the program are the material constants 
and relevant geometric quantities expressed in non-dimensional forms, 

namely, h, h c , l Q , & c , Vi . V, and V* • The solution procedure is as 

follows: Choose some initial value of w, and compute the roots of = 

y 2 t and also find all of the X p by solving Equations (5.38) and 
(5.39). Construct the 16x16 complex D matrix and check to see whether 


Equation (5.55) is satisfied. If not increment u and repeat the 
entire process till a zero crossing for the value of the determinant D 
is obtained. 

The solution for the undamped system (i.e., D R (w R ) = 0) is a good 

starting point for numerical Iteration. This solution can later be 
used as an approximate solution for, D R (w R +iO) + i (t/ R + 1 0) = 0 , 

and iteration is carried out by a bracketing method to arrive at the 
complex eigenvalue. 

Numerical Results and Discussion 

' The length and thickness of the graphite epoxy beams chosen here 
for study were: Z. = 12.25 cm and h = 0.25 cm. The storage modulus of 

the beam material as supplied by the manufacturer was 124 Gpa. The 
material loss factor of the beam was taken as = 0.004, : which was 

assumed to be a constant for all the modes. This corresponds to a 
damping ratio of 0.2% and is believed to represent closely the 
material damping of the graphite epoxy material as found from various 
experiments described previously In Chapter 3. The thickness of the 
adhesive layer was taken as 0.5 cm. The complex modulus of the 
adhesive material, which is epoxy resin in the present case was 
assumed to be 4(l+i0.04) Gpa which has a real part of 4 Gpa as 
supplied by the manufacturer. The complex shear modulus of the 
adhesive is assumed to be 1.4(l+i0.04) Gpa. Here tj 2 = i ] 3 = 0.04, 



which corresponds to a constant damping ratio of 2% for the epoxy 
material . 

The above quantities expressed In non-dimensional form were used 
In the numerical computation of the eigenvalues and eigenvectors. The 
ratio of the length of overlap to the beam length is called the 
overlap ratio. The Numerical results were obtained for overlap ratios 
of 0.2, 0.4, 0.6, and 0.8. Figure 42 shows the variation of natural 
frequency with overlap ratio for the first three modes of vibration. 
From this plot, it is clear that the natural frequency of the system 
Increases with an increase in the overlap ratio for all modes. Figure 
43 is a similar plot showing the variation of the system damping 
(total loss factor of the system) with overlap ratio. For the case of 

overlap ratio = 0.2, the non-dimensional frequencies u are: 2.36, 
13.12, and 27.24 for the first three modes. The corresponding system 
damping ratios are ^ , = 0.534%, £" a = 0.120%, and £, = 0.054%. It 

should be noted that the above values were obtained by assuming 
constant values of 0.2% and 2.0% for the material damping of the beam 
and the adhesive respectively. It is seen that, for this case, the 
system damping is higher than the material damping of the beam only 
for the first mode. For the other two modes, the system damping ratio 
appears to be lower than the material damping of the beam. This 
observation, however, is true only for the above case and cannot be 
generalized. To substantiate this point, let us consider a different 
case in which the damping ratio of the adhesive was assumed to be 5%, 



and the material damping of the beam was kept the same. The system 
damping ratios computed for this case are: ^ = 1.66%, * 0.33%, and 

£-3 = 0.16%. There were no significant changes in the natural 

frequencies of the system. As can be seen, the system damping ratios 
In this case are totally different from the previous case considered. 
Hence, In order to predict the system damping ratios, using the 
present model, accurate knowledge about the material damping ratios 
(of the beam and adhesive) and their dependence on frequency Is 
mandatory. 

Figures 44 through 46 show the predicted mode shapes of the first 
three bending modes for different overlap ratios. These were obtained 
by solving the homogeneous algebraic Equations (5.53) for each 
eigenvalue. 







Overlap Ratio 

Figure 43. Plot of Damping Ratio vs. Overlap Ratio 








Comparison with Experimental Results 


Table 5 shows both theoretical and experimental results of two 
graphite epoxy lap joint systems with overlap ratios of 0.2 and 0.5. 
All of the theoretical results were generated by considering simply- 
supported boundary conditions at the ends. 


Table 5 

Comparison of Theoretical and Experimental Results of the 
Graphite Epoxy Lap Joint Systems 


Overlap 

Ratio 

Mode No. 

Natural Frequency (Hz) 

Damping Ratio(%) 

Theory 

1 Expt* 

%Error 

Theory* 

Expt! 

Expt? 


1 

195.9 

207.0 

5 

0.53 

1.51 

0.20 

0.2 

2 

898.5 

986.2 

9 

0.12 

1.63 

0.16 


3 

1864.5 

1984.0 

6 

0.05 

2.29 

0.10 


1 

263.5 

270.6 

3 

0.13 

0.68 

0.43 

0.4 

2 

1148.4 

1075.5 

-6 

0.15 

1.43 

0.12 


3 

2492.7 

2684.4 

7 

0.20 

0.83 

0.15 


*Using Simply-Supported boundary conditions 
2 

Using Free Free boundary conditions. 


The specimens were prepared by bonding two similar graphite epoxy 
beams over the desired length of overlap using an epoxy adhesive. The 
epoxy resin was procured from CIBA-GEIGY Corporation. Much care was 
taken to obtain a good bond by properly curing the joint system in an 
oven. The dimensions of the unbonded beams were: length = 12.25 cm, 
width = 2.8 cm, and thickness = 0.25 cm. These dimensions and the 
material constants (Young's Modulus and Shear Modulus) were input to 
the computer program to predict the theoretical natural frequencies, 
modal damping ratios and mode shapes. 

The supports to simulate simply-supported boundaries at the ends 
were specially fabricated in the form of two separate triangular 
blocks with a knife edge on each, to support the beam, one from the 
top, and the other from the bottom. The two blocks were then clamped 
to a rigid base. In order to evaluate the experimental simulation of 
these supports, a trial test was first conducted on an aluminum beam. 
Excellent agreement was obtained between the theoretical natural 
frequencies (computed using simply-supported boundary conditions) 
andmeasured values. In fact, the percentage difference between the 
two results never exceeded 2% for the first four modes. Having 
established the validity of the above text fixtures, experiments were 
later conducted on graphite epoxy lap-jointed beams. An impact hammer 
with an attached force transducer was used for exciting the specimen 
and the response was measured using a mini -accel erometer (Bruel & 
Kjaer 4375). The frequency response (ratio of acceleration to force 
signals) was immediately computed and recorded on an FFT analyzer 



137 


(Bruel & Kjaer 2032). The modal parameters were then computed using 
the Structural Measurement Systems (SMS) modal analysis software. 

Referring to Table 5 it is seen that there is good agreement 
between the predicted values of natural frequencies and experimental 
data. The percentage difference between the two results is in the 
range of 3 to 9%. The measured damping ratio values in column 7 
(Table 5) are larger than the predicted values. This is presumably 
due to additional damping contributions at the end supports in the 
experimental data. In order to substantiate this, damping 
measurements were made on the same samples without the end fixtures. 
Column 8 of Table 5 shows these results. In this case, testing was 
done by simply mounting the sample directly on a shaker using a thin 
layer of wax. The resonance frequencies In this case were, of course, 
somewhat different from those obtained previously using the simply- 
supported boundary conditions. It Is Interesting to notice, however, 
that the predicted damping values in column 6 and the measure'd data in 
column 8 are of the same order. Next, theoretical mode shape is 
compared with experimental mode shape In Figures 47 for mode 1. Note 
that no attempt has been made to normalize the different amplitudes in 
the theoretical and experimental results. There is excellent 
agreement between the two mode shapes for the first mode as seen in 
Figure 47. The small discrepancies In the two results can be 
attributed to several assumptions made in the theoretical analysis. 
The major assumption is the use of a constant value with frequency of 
vibration for the modulus of the beam material. Also, in the 


V/ 



theoretical model it was assumed that the adhesive Is flexible and 
that the damping is caused by both shering and transverse motions of 
the adhesive. In practice, however, it Is very difficult to fabricate 
specimens to satisfy exactly the above requirements. 



Theory 



Experiment 



Figure 47. Comparison of Theoretical and Experimental Mode Shapes for 
Overlap Ratio = 0.2, Mode 1 



Experimental Study of Damping of Bonded 
Double-Butt Joint Specimen^ 


This section contains the results of experiments conducted on 
beam and jointed specimens in a vacuum chamber. Damping measurements 
were made on the following ten graphite epoxy samples which were 
procured from Sikorski Aircraft Company: 

i) Five identical graphite epoxy beam specimens with a double-butt 
joint in the middle as shown in Figure 48 and ii) Five identical 
graphite epoxy beam specimens of dimensions as in 1 but with no 
joints. 



Figure 48. Double-Butt-Jointed Specimen 


In each of the above experiments, the specimen was mounted at its 
center directly on the shaker using a thin layer of wax. This closely 
simulated free-free boundary conditions at the edges. Only odd 
numbered modes were excited, since the specimen was mounted at the 
center, which happens to be a node point for all even numbered modes 
for free-free boundary conditions. All the damping measurements were 
made using the usual improved half-power points method. 

First, measurements were made on the doubl e-butt- joi nted (DBJ) 
specimens. In each case, the specimen with the shaker was kept inside 
the vacuum chamber and damping ratios and resonance frequencies were 
measured for the first and third mode of vibration, before the air was 
pumped out. Then the measurements were repeated in a vacuum when the 
pressure Inside the chamber was about 0.1 mm of Hg. The test 
temperature remained ambient in both cases. Table 6 shows a summary 
of results for the five DBJ specimens tested under vacuum conditions. 
No significant change in the values of the damping ratio was noticed 
between the experiments conducted under normal atmospheric pressure 
and in a vacuum. This difference in the values of the damping ratio 
was in fact in the range of 0.01-0.02% during most of the trials. 

From Table 6, it is seen that the average value for the damping 
ratio of the DBJ specimen for the first mode is 0.114% and that for 
the third mode is 0.155%. The small discrepancies in the values of 
the resonance frequencies in the five samples could be attributed to 
a) slight differences in their dimension and/or b) slight differences 
in the exact location of the excitation point during mounting. 


The next set of experiments was conducted on the five beam 
specimens with no joints (having the same dimensions as the DBJ 
samples) under identical environmental conditions as maintained with 
the DBJ specimens. The results are tabulated in Table 7. Here, it is 
seen that the beam has an average value for the damping ratio of 
0.107% for the first mode and 0.138% for the third mode. 

By comparing the results of Table 6 and 7, it is seen that, 
although the damping ratio values for the jointed specimens appear 
somewhat higher than those of the beam specimens, this difference is 
very small. This observation is also true for the previous damping 
results of lap-jointed specimens with free-free boundary conditions. 

Much care was exercised in maintaining Identical environmental 
and other conditions in all of the above experiments. The damping 
measurements were made precisely by zooming in on each mode thus 
increasing the frequency resolution of the analyzer and minimizing the 
influence of external noise. Figure 49 is a plot of the, frequency 
response data for a DBJ specimen for baseband excitation In the 
frequency range 0-1.6 kHz. The corresponding zoom frequency response 
plots for modes 1 and 3 are shown in Figures 50 and 51 respectively. 
Figures 52 through 54 show similar baseband and zoom frequency 
response plots for a beam sample. 



Experimental Results of Double-Butt-Jointed Specimens 


Sample No. 

Mode 1 

Mode 

3 


f <Hz) 

Ci,M 

f n (Hz) 


1 

258.2 

0.11 

1371.2 

0.14 

2 

257.0 

0.10 

1375.9 

0.13 

3 

260.0 

0.13 

1381.7 

0.17 

4 

256.3 

0.10 

1368.5 

0.17 

5 

261.3 

0.12 

1380.1 

0.16 

Average 

258.6 

0.11 

1375.5 

0.15 



Table 7 



Experimental 

Results of 

Beam Specimens with 

No Joints 

Sample No. 

Mode 1 

Mode 

3 




f n (H z > 


1 

341.1 

0.11 

18C0.0 

0.15 

2 

337.3 

0.13 

1782.2 

0.15 

3 

340.5 

0.11 

1794.3 

0.12 

4 

330.4 

0.11 

1749.2 

0.12 

5 

334.5 

0.10 

1769.0 

0.15 

Average 

336.8 

0.11 

1778.9 

0.14 
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Figure 49. Frequency Response Plot of the Double-Butt-Jointed 
Specimen Using Baseband Analysis 



Figure 50. Frequency Response Plot of the Double-Butt-Jointed 
Specimen Using Zoom Analysis for Mode 1 



Figure 51. Frequency Response Plot of the Double-Butt-Jointed 
Specimen Using Zoom Analysis for Mode 3 
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Fiaure 52. Frequency Response Plot of the Been, Specimen Using 
Baseband Analysis 



Figure 54. Frequency Response Plot of the Beam Specimen Using 
3 ~ * Zoom Analysis for Mode 3 
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VI. SUMMARY AND CONCLUSIONS 

In this project an attempt has been made to evaluate the damping 
capacity of graphite epoxy materials and structural joints. The 
proj ect was sponsored by the NASA George C. Marshall Space Flight 
Center, Huntsville, Alabama under contract No. NAS8-36146. The work 
^involved a systematic investigation of tjie damping ratio of different 
composite specimens and bonded joints^under normal atmospheric 
conditions and in a vacuum environment. Free and forced vibration test 
methods were employed for measuring the damping ratios. The effect of 
edge support conditions on the damping value of a composite tube 
specimen was studied by using a series of experiments performed on the 
specimen with different edge supports. It was finally : found that 
simulating a free-free boundary conditions by having no constraints at 
the ends gives the lowest value of the material damping of the 
composite. The accuracy of the estimation of the damping ratio value 
was Improved by using a curve-fitting technique on the response data 
obtained through measurement. 

The effect of outgassing (moisture desorption) on the damping 
capacity was determined by measuring the damping ratio of the tube 
specimen in a vacuum environment before and after outgassing had 
occurred. The effects of high and low temperatures on the damping was 
also investigated by using a series of experiments on tube and beam 
specimens. Special experimental set-ups were developed to accomplish 



—this objective^ An analytical model to study the vibrations of a 
bonded lap joint system was formulated. Numerical results were 
generated for different overlap ratios of the system. These were 
compared with experimental results. 

In order to determine the Influence of bonded joints on the 
material damping capacity, experiments were conducted on bonded lap- 
jointed and double-butt-jointed specimens. These experimental results 
were compared with simple beam specimens with no joints.- _ v "‘ 

The following conclusions may be drawn from the results of this 
investigation. 

1) The damping ratio value is strongly dependent on the end support 
conditions used to support the specimens. It was observed that 
the end fixtures usually offer a large contribution to the damping 
capacity of the material. 

2) The material damping ratio value of the composite tube is in the 
range of 0.1% - 0.15% for a frequency range of 500 Hz to 560 Hz 
under normal atmospheric conditions. The first bending mode of 
the specimen occurs in this frequency range. The damping ratio 
value seems to increase with higher modes and has a value of about 
0.36% for a modal frequency of about 2840 Hz. 

3) From the experiments conducted on the tube specimen in a vacuum, 
it may be concluded that the contribution of air damping to the 
material damping is between about 0.03% and 0.04%. The material 
damping ratio of the tube specimen in a vacuum is thus in the 
range of 0.06% and 0.11%. The effect of outgassing or moisture 
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desorption on the damping ratio value of the specimen Is not very 
significant. The change in damping ratio value due to outgassing 
of the tube specimen for a period of about one month was observed 

V 

to be 0.016%. The percentage change in the weight of the tube 
from its original weight was 0.9%. No significant change in the 
damping ratio with the vacuum pressure inside the chamber was 
observed. 

4) It is observed that changes in temperature have a significant 
influence on the damping capacity of the composite material. From 

V 

. the experiments conducted on tube and beam specimens, It is 

observed that the damping of the material increases with 

temperatures above the atmospheric temperature and has a value of 
w 

about 3% near the material glass transition temperature. There 
is, however, a small decreasing trend in the values of the damping 
with temperatures below the normal atmospheric temperature. This 
is In agreement with the data available in the literature on other 
types of composite materials. 

5) The analytical model described in Chapter 4 can be used to predict 
the natural frequencies and the modes shapes of a bonded lap joint 
system for free vibration. The model can also be used to predict 
the system modal damping values by properly choosing the material 
damping values of the beam and the adhesive. Good agreement 

1 * 

between numerical and experimental results for the modal 

frequencies and mode shapes was obtained. The percentage 
v 


di fferrence between the values of the theoretical and experimental 
natural frequencies Is In the range of 3 to 9%. 

6) From the numerical and experimental results obtained on lap- 
jointed beams, it is clear that the natural frequencies of the 
system Increase with an Increase in the length of overlap* 

7) The increase In the value of the damping ratio due to the presence 
of bonded joints in the system does not appear to be very 



VII. RECOMMENDATIONS FOR FUTURE WORK 


The material damping of the graphite epoxy specimens determined 
so far appears to be very low. This low value of damping may result 
in a significant truss system displacement amplitude of the telescope 
truss system, which, in turn may lead to a blurring of the optical 
system image. The total truss system damping, of course, will be be 
slightly higher because of the effect of the joints. But the precise 
effects of the bonded joints on the system damping is not known. 
Unbonded joints are known to provide significant contributions to the 
damping. More work is thus needed in this area to determine if the 
bonded joints can be relied on to increase the value of damping. This 
work should include both analytical and experimental work on other 
types of joints as encountered in space structures. Some of these 
joints are tubular joints, bolted joints etc. Typical configurations 
of tubular joints are shown in Figure 55. 

Another area which needs attention is the dependence of damping 
with frequency. The material damping at very low frequencies must be 
determined since many large space structures are found to possess very 
low fundamental resonance frequencies which are generally less than 
about 10 Hz. The dependence of damping with displacement amplitude is 
also important when the system vibrates at such low frequencies. 

The study on the effect of moisture content on the damping 
capacity of these materials is not complete in many aspects. The 


moisture absorption and desorption capacity of these materials Is not 
well known. The percentage moisture content of the material with 
respect to time, and the exact glass transition temperature need to be 
determined. The moisture absorption capacity of these materials is 
expected to be different if the material is soaked in water or exposed 
to humid air. These two approaches are expected to produce entirely 
d "j -f f e p q n t results. It is thus important to examine both the 
approaches in the future to determine the exact influence of moisture 
absorption and desorption on the damping capacity of these materials. 

Work is also needed in the prediction of damping of the graphite 
epoxy material. Some work has been done by other investigators in 
predicting the damping of laminated composite plates using a finite 
element approach. There is a considerable scope for research in this 
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ILLUSTRATION OF THE CURVE-FITTING TECHNIQUE 
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ILLUSTRATION OF THE CURVE-FITTING TECHNIQUE 
The half-power points method of estimation of the damping ratio is 
not always accurate. It relies on how well one can read the frequency 
values for the half-power points from the frequency response data. It 
is seldom possible to estimate the half-power points without making any 
interpolation in the frequency values. The curve-fitting technique 
described in Chapter 2 can be used to overcome the above problem with 
the half-power points method. In the curve-fitting method, the entire 
data near a resonance are used for the estimation of the damping ratio, 
as opposed to three points in the half-power points method. Thus, a 
more accurate estimation of the damping ratio can be obtained. This 
Method has been found to be especially well-suited and more accurate 
than the half-power points method for the case in which the frequency 

response curve is contaminated with noise. 

In the present illustration, the damping ratio of the tube specimen 
was estimated by both the methods and the results are compared. In the 
curve-fitting technique, all three quantities, namely, the magnitude, 
'the real part, and the imaginary part of the receptance function were 
utilized separately to estimate the damping ratio value. Figure 56 
shows a comparison of the plots for the magnitude of the frequency 
response function obtained from both the measurement and the curve- 
fitting algorithm. Figures 57 and 58 show similar plots for the real 
part and the imaginary part of the frequency response function. All 



these data are for the tube specimen with a half ring as the connecting 
element. It may be recalled that the damping ratio value for the data 
using the half-power points method is 0.130% and the resonance 
frequency is 552.50 Hz. Now, using the curve-fitted response the 
damping ratio is estimated to be 0.128% and the resonance frequency to 
be 551.45 Hz. From the real part of the receptance curve using the 
curve-fitting method the damping ratio value is estimated to be 0.1294 
and the resonance frequency to be 551.36 Hz. From the imaginary part or 
the frequency response function and the curve-fitting method the 
damping ratio value is found to be 0.121% and the resonance frequency 
is 551.63Hz. All of the above results are summarized in Table 8. It 
Is clear that there is excellent agreement In the results from the 
different approaches used. 

In the present illustration the measured frequency response curve 
is close to an 'ideal 1 response plot of a single degree of freedom 
system. It shows a well-defined peak and the frequency range 'chosen is 
small. It was thus possible to easily locate the resonance frequencies 
and the half-power points. Hence, the damping ratio estimated from the 
half-power points agrees well with that of the curve-fitting method. 
But, in cases where the frequency response curve is not so 'ideal 1 as 
in the present case, it is always advantageous to estimate the damping 
ratio value using the curve-fitting method than the half-power points 


method . 



Table 8 
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Comparison of Damping Data from Different Estimation Methods 




Curve-fitting data 

Measured data 


From magnitude 
plot of 
receptance 

From real part 
plot of 
receptance 

From imaginary 
part plot of 
receptance 

From half- 
power points 
method 

Damping 

ratio 

0.128% 

0.129% 

0.121% 

0.130% 

Resonance 

frequency 

(Hz) 

551.45 

551.36 

551.63 

551.50 


; .V 




V 
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Receptance- Real Part Receptance- Magnitude m/)l 


V 





Figure 56. Comparison of the Curve-Fit and Experimental Data of the 
Magnitude of the Receptance 



FREQUENCY-Hz 


1. Measured data 

2. Curve-fit data 


Figure 57. Comparison of the Curve-Fit and Experimental Data of the 
Real Part of the Receptance 





APPENDIX B 

EQUATION OF MOTION OF STANDARD LINEAR 
SOLID MODEL 





EQUATION OF MOTION OF STANDARD LINEAR SOLID MODEL 
The standard linear solid model Is shown in Figure59 (a). It 
consists of a series combination of an elastic spring Kj , and viscous 

dashpot c, connected in parallel with another elastic spring K. When 
Kj approaches infinity, the standard model reduces to the well~known 

Kel vin-Voight model as shown in Figure 59(b). In the other special 
case, when K is identically zero, the standard model reduces to the 
Maxwell model. The equation of motion for the standard linear model is 
developed as follows. 

Let Pi and P 2 be the forces in springs Kj and K, such that 


'• 


P = 

Pi + P 2 . 

(Bl) 

Referring to 

Figure 

59(a), Let 

y 3 be the displacement of 

the node 

connecting K 2 

and c. 

The forces F 

\ and P 2 are given by 




Pi - Ki 

(yi-y 3 ), and 

(B2) 



P 2 = K 

• 

« 

1 

^>1 

(B3) 


Differentiating Equation (B2) w.r.t. time t, we get, 

3 Pi 3Yi ay 3 

at - = Kl at at ' * 


(B4) 


The force in the dashpot is, 

Pi 


(3y 3 3y 2 ^ 

aF " aF 


(B5) 


From Equations (B4) and (B5) , 
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w 


f 


w 


3y x l 9 Pi Pi Q yj 
aF = iqatT + c" + at ' 

ay, 3y 2 1 ePj Pi 

1 * “ • F _ — — ■ • + * , 

at at K,at c 

From Equation (Bl) we have P » P, + P 2 • therefore, 

8P 9P X 3Pj 

at 

Hence Equation (B7) becomes, 

»y» 1 r ,_ 

at "at " Kj L ^at " at 


at at 

3P 3P a 1 


(B6) 

(B7) 


(B8) 


(B9) 


Substituting for P 2 and from Equation (B3) , we obtain, 

at 


3Yi ay 2 _ 1 SP v *Vx_ 

at at K t *• at 


- ^ - v [ ST - K(ar - er » + c [p - K(y -- y ’ )] - 


(BIO). 


Rearranging, 


1 3P P ay 

KjFt + c = fe 


8y 2 K K 

at ) [1 + k, ] + c 


(Bl 1) 


This is the equation for the general linear solid model. In the limit 
when K, approaches infinity, the standard linear solid model becomes 

the Kel vin-Voight model as shown in Figure 59 (b) . Taking the limit of 
Equation (Bll) as K x approaches infinity, we get 

By ay , 

p- c (ar--ar> + K(y ‘- y ’ ) - 


(B12) 
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(a) Standard Linear Model 



(b) Kelvin-Voight (c) Maxwell Model 

Model 


Figure 59. Visccelastic Models 
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VIII. IMPROVEMENTS IN DAMPING MEASUREMENTS 


In this chapter some improvements that have been considered in making 
accurate damping measurements are discussed. Since our main aim is to measure 
the damping of the graphite epoxy material, the damping from the supports of 
the specimen and other damping, such as that from the connecting cables and 
surrounding air are to be avoided. As discussed in the previous chapters, the 
best way to avoid damping from supports is to use a free-free support 
condition for the beam and excite it at its center through an impedance head 
mounted on a shaker. The input force signal to the beam and the output 
acceleration signal from it are fed through two charge amplifiers to the Dual 
Channel Signal Analyzer; Then the frequency response either inertance A/F, 
mobility V/F or receptance Y/F is obtained. 

Preliminary Check of the Measured Results 

The accuracy of the measured results may be evaluated from an observation 
of the frequency response curve and the measured resonance frequency by 
•comparing them with the values obtained by theoretical analysis. 
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The Shape of Frequency Response Curve 

A mathematical analysis is described herein which would help in 
predicting the frequency response curve for simple beams studied (the beams 
without joi nts) . 

For a simple beam with material and air damping as mentioned above, its 
displacement y can be obtained by superposition of displacements y i 

corresponding to its infinite normal modes of vibration X i (x) . 


y(x,t) = [ y-j = E ^ (t) x i ( x ) 

i=0 1=0 

where normalized displacement function [121, 122] 

X 0 (x) - — (rigid-body mode) 

VI 


and 


X, 


X, 


X i ( x ) = [(coshj^x + cos-jx) - o’ i (sinh— j-x + sin— yx)] 

coshX^ _ cosX. 

(i = 1,2,3 ...» sinhX.. - sinX^ 


( 8 . 1 ) 


( 8 . 2 ) 


(8.3) 


are the solutions for the differential equation (8.4) with boundary conditions 
as given in Eq. (8.5) . 

(8.4) 


tA ax« 


+ = 0 


?li| . Siii , . o 

3x 2 l x=0 3x 2 ! x=l 


£*i B 2l| = o 

3x 3 'x=0 3x 3 1 X=1 


(8.5) 



According to W.E. Baker et ai [123], for a beam with material and air 
damping and with small amplitude vibration,, its differential equation of 
motion for free vibration analysis is 


311 + g ft 91X_ + £ ^ + a 3H = 

a# + ^ EI ax^at 1 at p %v 


( 8 . 6 ) 


Where the meaning of /? can be found from Eq. (5.15) of this report. From Eqs. 
(8.1) and (8.3), X i (x) can be separated from Eq. (8.6) and the differential 
equation for the time function ^ (t) (irO) is obtained. 


+ 2a i + w* <j>\ = 0 (irO) 


where 


2 /XV El 

w i = pA 


2a i = V E + JK) 


(8.7) 


( 8 . 8 ) 


(8.9) 


By using the principle that the total virtual work should be zero, an equation 
similar to Eq. (8.7) can be written for a free-free beam with damping, 
supporting and exciting conditions as mentioned above. 

jj. + 2 a.fj + w* ^ X- (^)coswt + glX^ (|) (i/0) 1 (8.10) 
where F 0 coswt is the exciting force and g is the acceleration due to gravity. 


However, for the rigid-body mode, the differential equation is 


! y 0 3y 0 

— = F 0 coswt - C 


( 8 . 11 ) 


because there is no elastic force applied on the beam. From Eq. (8.11) the 
differential equation for is written as 


S 0 + &<f> = coswt = coswt, 

/;A1X 0 An/1 


( 8 . 12 ) 
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where 


w 







a = 


(8.13) 


/iAl 

According to S. Timoshenko et al [122], the solution with the initial 
conditions 

y| t -o = 0 

(8.14) 

= 0 

9t 1 t=0 


for Eqs. (8.12) and (8.10) are respectively as follows 


-at 


, ,.s _ 0 r cosjVt+j} e i 

" pM] L‘ „ 2 ..2, .* , JJ' 


(8.15) 


uV(u +a ) ui +a 


#,(*> - %1X,(J)[1 - £ e’ 0 !* cos(« d t- 7l .)] + 


F.X j (i)cos(ut-« j ) 


2 2,2 


where 


) +(2a^u) ] 


FoXj^e'V 


[(w*-w*)cosw d t + (^ + ^ )^sinw d t] (lVO) (8.16) 


«. *, G i_, 


^A[ («*-«*)* + (2^6/)*] 


}=tan _1 (|), 


(8.17) 


and 


7 .«tan _1 (j^) (iVO), 
d 


. 2a, w 

5.= tan -i ( i T ) (1*0), 
! w. -u 


d = V(w 2 -a 2 ) (irO) 


(8.18) 
(8. IS) 


( 8 . 20 ) 





Since X. (^)= 0 (i=2,4.6, ...«) » for the center of the beam, the displacement 
and the steady-state displacement of the above beam are as follows 


00 _ 

y(j.t) = ?J 0 (t)x o (l) + l (t)x i (^) 


( 8 . 21 ) 


1=1,3, 5 


and 


fl 


, 2 /I 


y(*.t) = - 


cos (ut +$} 


l , , c { 9 5T 1X 1 ( 2 ) + 


, ] F 0 X i (^)cos(«t-fl i ) 


2 „ 2 


} ( 8 . 22 ) 


^ Al wv'(u 2 +a 2 ) T=1,3,5V. 1 c pM[{u^- u ) + (2 a^u) ]' 

Where the first term is the displacement for the rigid-body mode due to the 
exciting force, and after the summation sign the first term is the static 
displacement due to gravity and support, and the last term is the steady-state 
displacement due to the exciting force. 

The corresponding velocity and acceleration at the center of the beam are 
written as 


9 v fl tl - — r s 1n(tft-h0. _ r 

at y{ 2 ' Z) ~ pA] t„.«, V 


IX 2 (^)wsin(wt-0..) 
2 2““2 


V(w + a ) i=l,3,5... >/[(w.-w ) +(2 a^w) ] 


,} 


(8.23) 


0 2 


5y 2 


. F 0 ui ' r 1 V 

y(£,t) = ttt { — 2 — rr cos(ut+f) - [ =1 3 5 t 

2 P A{ V(w + a ) i 4(w,- 


IX* (^)w a cos(wt-0.) 


'/[(w.-w) + (2a- u) ] 


} (8.24) 


By expressing Eqs. (8.23) and (8.24) in the complex form we have 






-l 


1X 2 (1) W e“ J ®1 


^ A1 V(w‘+a 2 ) i=l,3,5...7[(w 2 V) 2 + (2a 1 w)“] 

F 0 jt I) 00 IX. (i)w e ^i 

../l 0 rwe J r r i'2' 1 

yU.W ~ /,Ai \ 2 2 " L r 2 2 2 ,2-|j 

V(w +a ) i=l,3,5... >/[(w. - u ) + (2a^ui) ] 


8 2 


ay 


} e^ wt "? (8.25) 


Jwt 


(8.26) 


or 
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/\ 

a , F 0 a a j (wt - f + tan" 1 

| 4 ^MT V(A “ + E « )e 


A, 


F 0 W 2 2 


= Tfli ^(A 0 + B 0 )e 


_i ° 

j (wt + tan g-) 


at 


/>A1 


where 


« _ cos l 

n 0 .22 


l 


!X^.(^)w cos0 i 


V(w +a ) i =1 ,3,5... '/[(w.. - u ) + (2 a^w) ] 


co 

B„ . 51M-5- + l 


1 X* (^)w sin#.. 


V(w Z +a 2 ) i=l, 3, 5. . •v'[(w^ -u ) + (2a. u) ] 


(8.27) 


(8.28) 


(8.29) 


(8.30) 


or 


2 ] 2 2 

CO IX. (i)w (w. - hi ) 

a _ ^ r 1_£ J 

A o - 2 2 L 25.2 ! 

u +a 1— 1,3,5... (w.-w ) + (2a. w) 


Bo » 4-r + l 


lX*(l)w(2a.tf) 
2 272 


(8.31) 


(8.32) 


u +a 1=1, 3, 5... (wj - w ) + ( 20 ^) 

Then the magnitude of the frequency response, the mobility and inertance are 
respectively as follows: 


||h = ll- — =h 7(a « + B ’ J 

r° *° 3 1 ^ 

and 

Ie I = [f-^ry(j.t)l “h v(a " + b °> 

r 0 r 0 Qt ^ 

where 

M=/?A1 


(8.33) 


(8.34) 


(8.35) 



Figs. 60 and 61 respectively show the theoretical frequency response 
(mobility and inertance) of a graphite epoxy beam with the following 
parameters . 

1 = 2.54xl0 -1 m (length of beam) 
b = 2.84xl0 _2 m (width of beam) 

r> 

h = 2.65x10 m (thickness of beam) 

E = 1.2X10 11 Pa (storage modulus of elasticity for graphite epoxy) 

3 3 

p = 1.6x10 kg/m (density of graphite epoxy) 
q = 0.001 (damping ratio of beam) 
a = Wif 
a i = 070) 

When plotting these frequency response curves from Eqs. (8.33) and (8.34), 
only four leading odd modes of vibration were used in Eqs. (8.31) and (8.32). 

By comparing the experimental frequency response of the mobility and 
Inertance shown in Figs. 62 and 63 with the theoretical ones shown in Figs. 60 
and 61 respectively, it is possible to estimate the accuracy of the measured 
results. Otherwise, if there are significant variations between the two 
curves, it is better to check the problem before the next measurements are 










Figure 62. Experimental Frequency Response of Nobility of the Beam and 
Skeleton Check 
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Figure 63. Experimental Frequency Response of Inertance of the Beam 



The Skeleton of Mobility Plot [25] 


For a structure more complicated than a simple beam, it is difficult to 
predict its frequency response. Therefore, the utility of observing the shape 
of the frequency response curve is limited. Another technique used for 
checking the measured results is called skeleton check and is shown in Fig. 
w 62. It is suitable for any structure with free-free support condition, and it 

should be used for the frequency response curves which are plotted in 
logarithmic scales. If there is a tendency that the ±6dB/0ct straight lines 
V go down with an increase in frequency, the measured results may be all right. 

We can find our results to agree very well with this method of checking. 

v The Accuracy of the Resonance Frequency 

For a lightly damped simple beam it is easy to predict-its natural 
•w frequency. Hence, a good agreement of the measured resonance frequencies with 

the predicted natural ■ frequencies also indicates a good measurement. 
Obviously if the measured resonance frequencies are found to be lower than the 
— calculated ones there must be a load effecting the beam vibration. 



Improving the Mounting of Beams 


It is found that the mounting of the beam specimens under test play an 
important role in damping measurements. Hence time and effort has been spent 
for selecting the best mounting method. The advantage of supporting the beam 
at its center is that it prevents excessive damping from the supports as in 
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the case of a cantilevered or a simply supported beam. Unfortunately, the 
di sadvantage ’for a center supported beam is that the center point is the 
antinode of the odd mode vibration of the beam, and that a load at the 
antinode will produce a notable impact. If the load is a mass, it will only 
lower the resonance frequencies; if it is a resistance, it will not only lower 
the resonance frequencies but also increase the damping [124]. 

After extensive trials with wax, stud, glue and tape, the double sided 
adhesive thin tape has been chosen for mounting the specimens. Wax is not 
ideal, because a mounting with less wax tends to be too loose while the 
specimen is vibrating, and if too much wax is used, the measured damping ratio 
will be very high. Although the stud mounting is stable enough, it tends to 
decrease the resonance frequencies and increase the damping ratio value 
significantly. Glue mounting is also stable, but the surface area of contact 
between the beam and the impedance head is often too large, hence the measured 
damping ratio value will tend to be higher than that of the specimen under 
test. 

In order to get a rather stable and minimally damped support- for making 
damping measurements, a narrow piece of carpet tape about 1 mm in width was 
stuck at the center of the beam perpendicular to the length of the beam. 
Since the tape is very thin, it was necessary to keep the beam balanced. This 
was done by ensuring that the gap between the beam and the impedance head on 
either side of the tape mounting is the same. 

Selection of Excitation 


There are several kinds of signals generated by the B&K Dual Channel 
Signal Analyzer Type 2032. Pseudo-random noise is a self-windowing function, 
and it is often used as the excitation in measuring the frequency response of 



fairly linear systems, particularly for lightly damped resonances, because the 
self-windowing functions generate no leakage while obtaining the Fast Fourier 
Transform. Without a self-windowing excitation, even with a Hanning Window, 
energy could leak from a lightly damped resonance into adjacent spectrum 
lines. The resulting spectrum would show greater damping than it actually 
exists. Furthermore, pseudo-random noise gives the frequency response much 
faster when compared to the random noise excitation. Because the latter is a 
true noise, it must be averaged for several time records before an accurate 
frequency response can be determined [125]. 

Another aspect of the excitation is its amplitude. The rule of thumb for 
selecting a suitable amplitude is that the excitation must be strong enough to 
get the needed signal to noise ratio, and must not be so strong as to damage 
the stability of mounting and linearity of the system. 

Elimination of Noise in Measurements 
Electrical Noise 

Electrical noise may be eliminated by properly grounding the measuring 
instruments according to their specifications. It is necessary to point out 
that the housing of the impedance head should be isolated from the grounded 
shaker with a specially made electrically isolated stud and a mica washer. 

There is often an electrical noise at the center frequency of the 
frequency span in the frequency response measurements made by B & K 2032 
Analyzer. In order to avoid this type of noise, it is suggested that the 
resonance frequency of the specimen be not selected as the center frequency of 


zoom measurements. 



Mechanical Noise 


Sometimes vibration coming from the foundation of the shaker very much 
effects the measurements. A piece of soft foam may be inserted beneath the 
shaker. This acts as a cushion and was found to be quite helpful in 
preventing such disturbance. But in the case of measurements taken inside a 
vacuum chamber, since the foam is often full of dust which is harmful for i.he 
vacuum pump, a vibration isolater with very low resonance frequency (just a 
few Hertz) was used in place of soft foam. 

Another problem is the force or vibration coming from the cables 
connected to the impedance head. Compared with the weight of the specimen, 
the disturbance of the cables is rather heavy, although the cables are soft 
and light. Keeping the cables as straight as possible to prevent them from 
bending and putting the cables on a cushion to damp out their vibration may 
solve this problem. But when working inside the vacuum chamber, it was found 
difficult to do so, because the space in the bell jar was very.-restrl cted. 
Therefore, it is suggested that softer and lighter cables be used while making 
measurements in a vacuum chamber. 

Further Check of the Measured Results 

After the above problems having been considered, the accuracy of the 
measured results improved to a great extent. In order to make sure that the 
experimental results were reliable, further check was always conducted during 
the measurement procedure. The signal to noise ratio and coherence between 
the two channels were observed for each measurement. The optimum requirements 
were found to be a signal to noise ratio of over 20dB, a coherence of over 




0.91, and a very small difference between the frequency responses H2 and HI 

measured by the B&K Analyzer. At the same time, the repeatability of several 
successive results was also observed. It was found that as the signal to 
noise ratio, coherence and the value of the damping ratio decreased, the 
difference between H2 and HI and the resonance frequency increased quickly. 
The reason was most often found to be related to the mounting of the beam 
which probably got slack. 


Experimental 


of Jointed 


)ecimens 


This section contains the results of experiments conducted on beam type 
specimens with and without joints. The joints studied herein were made from 
pieces of graphite epoxy beam type specimens. Damping measurements were made 
on the following type of jointed specimens. 1. Bolted Single Lap Joint (Fig. 
64), 2. Bonded Double Lap Joint (Fig. 65), 3. Bolted Double Lap Joint (Fig. 
66), 4. Bonded-Bol ted Double Lap Joint (Fig. 66), 5. Bonded Scarf Joint 
(Fig. 67), and 6. Bolted Butt Joint (Fig. 68). The measured damping ratios 
for these jointed specimens are shown in Table 9. The damping ratio value of 
specimens without joints was also measured under identical conditions when 
measurements were made on each of the jointed specimens. All the values of 
damping ratio reported in Table 9 were obtained after conducting about ten 
trial runs on each specimen. This was done to ensure repeatability of the 
results. Also each time a measurement was made it was checked for accuracy 
using the methods described earlier. 

The following conclusions may be drawn from the summary of results 
presented in Table 9. 
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The damping ratio of a graphite epoxy beam without joints is of the order 
of 0.05% to 0.06%. Compared to the earlier results listed in Tables 3,4 and 
7, the latest data are about one half of the previous ones. According to 
Plunkett [128], if measurements are taken on similar specimens, at the same 
stress level and by comparable methods, usually the lowest reported value of 
damping is most likely to be closer to the correct value. The actual material 
damping will probably be lower than any of the values that are measured. 
Also, there might have been some little differences in the material used for 
making the specimens whose damping values are reported in Tables 3,4 and 7 and 
those values reported in Table 9. 

The damping ratio value of jointed beams is found to be higher than that 
of the specimens without joints. This conclusion also differs from the 
earlier one described in pages 142 and 149. 

It appears that the damping ratio values tend to increase when the beams 
are bolted together than when they are just bonded together. The damping 
ratio value for a double lap joint is higher than that of a single lap joint. 

Lastly the difference in the damping ratio values measured under 
atmospheric conditions and in a vacuum (0.1mm Hg) seems to be not very 
significant. The possible reason is that the damping measurements might have 
been affected by the connecting cables used for the impedance head. 




i th Bolted Single Lap'Joints 














Damping of Various Graphite Epoxy Beams 



Beam 

Mode 

Resonance 

Frequency 

Damping Ratio 

No. 

Type & Size 

No. 

Atmospheri c 

0.1 mm Hg 

Atmospheric 0.1 mm Hg 


Without Joints 

1 

352Hz 


0.06% 

1 

254x28. 4x2. 65mm 


356Hz 

357Hz 

0.05% 0.04% 



3 

1871Hz 


0.04% 


Without Joints 

1 

472Hz 


0.06% 

2 

200x15x2. 5mm 






Bolted Single 

1 

343Hz 


0.11% 


Lap Joint 


344Hz 

344Hz 

0.11% 0.10% 

3 

254x28. 4mm 






Fig. 64 

3 

1618Hz 


0.12% 


Bonded Double 
Lap Joint 

4 190x15mm, 1 1028Hz 0.09% 

Fig. 65 


Bolted Double 
Lap Joint 

5 300x15mm, 1 1020Hz 0.11% 

Fig. 66 


Bonded-Bol ted 
Double Lap Joint 

6 300x15mm, 1 1021Hz 0.22% 

Fig. 66 


Bonded Scarf 
Joint 

7 200x15x2. 5mm 

Fig. 67 

1 

509Hz 

0.16% 

Bolted Butt 

1 

222 

0.51% 

Joi nt 




8 254x28. 4mm 

3 

1259 

0.11% 

Fig. 68 
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